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1 INTRODUCTION 



1.1 ftnrnnal MnHpls 



ABSTRACT 

We present spectral fits and timing analysis of Rossi X-ray Timing Explorer observa- 
tions of GX 339—4. These observations were carried out over a span of more than two 
years and encompassed both the soft/high and hard/low states. Two observations were 
simultaneous with Advanced Satellite for Cosmology and Astrophysics observations. 
Hysteresis in the soft state/hard state transition is observed. The hard state exhibits 
a possible anti-correlation between coronal compactness (i.e., spectral hardness) and 
the covering fraction of cold, reflecting material. The correlation between 'reflection 
fraction' and soft X-ray flux, however, appears to be more universal. Furthermore, low 
flux, hard state observations — taken over a decline into quiescence — show that the 
Fe line, independent of 'reflection fraction', remains broad and at a roughly constant 
equivalent width, counter to expectations from Advection Dominated Accretion Flow 
models. All power spectral densities (PSD) of the hard state X-ray lightcurves are 
describable as the sum of just a few broad, quasi-periodic features with frequencies 
that roughly scale as coronal compactness, £ c , to the —3/2 power. This is interpretable 
in a simple, toy model of an efficient spherical corona as variations of £ c oc i?t, where 
i?t is the 'transition radius' between the corona and an outer thin disc. Similar to ob- 
servations of Cyg X-l, time lags between soft and hard variability anti-correlate with 
coronal compactness, and peak shortly after the transition from the soft to hard state. 
A stronger correlation is seen between the time lags and the 'reflection fraction'. These 
latter facts might suggest that the time lags are associated with the known, spatially 
very extended, synchrotron emitting outflow. 

this type usually assumed that t he ADC sandwi ched a clas- 

sical Shakura fe Sunyaev (1975) accretion disc (Sunyaev & 

Triimper 197£ ; Haardt & Maraschi 1992 ) . It has generally 



been assumed that this coro na is produced via magneto- 



Most g alactic black hole candidates (BHCs) show transitions hydrodyna mical instabilities (Balbus & Hawley 1991 



between spectrally soft states and spectrally hard states. 
Hard spectral states appear to occur at bolometric lumi- 
nosities 5% of the Eddington limit, I/Ed d, while soft s pec- 
tral states occur at higher luminosities (Nowak 1995). A 
possible interpretation of this luminosity sequence is that 
state transitions are related to changes in the mass flow, M, 
through the accretion disc. The nature of the transition be- 
tween states, and, equally importantly, the geometry of the 
accretion flow within any given state, is still uncertain and 
a matter of vigorous debate (see Done 2001 for a review). 

Most current models favor Comptonization as the ma- 
jor physical process to produce the observed hard state X- 
ray spectrum (although see Markoff et al. 2001). In these 
models, soft photons from an accretion disc are Compton- 
upscattered in a hot (kT c ~ 150 keV), thermal electron 
plasma, the accretion disc corona (ADC). Early models of 



stone 



et al. 1996). It has been shown, however, that such a sand- 



wich geometry does not reproduce the observed X-ray spec- 
tra. The temperatures required to produce the hard power 
law cannot be reached since reprocessing of ~ 50% of the 
hard Comptonized radiatio n results in a very l arge Compton 
cooling rate of the corona (Dove et al. 1997a). 



In recent years, a large set of different ADC geome- 
tries have been proposed in order to circumvent the "heating 
problem" . All of these models have small covering factors for 
the ADC or attempt to otherwise reduce the fraction of re- 
processed hard radiation returning to the ADC as soft radia- 
tion. Models discus sed include: patchy, static, coronae above 



the accretion disc (Stern et al. 1995 



Gierliriski et al. 1997 



Poutanen & Fabian 1999); in- or outflowing coronae where 



the amount of reprocessed radiation is decreased due to rel- 
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Table 1. Observing log of the observations of GX 339—4. 



Obs. 


Date 


Exp. a 


3-9 keV 6 


9-20 keV 6 


20-200 keV c 


Optical 


Radio 




(y.m.d) 


(s) 




1n _ a —O —1 
10 erg cm s 




TV AT A 1 TV T 

Mb, My, Mj 


mjy 


P20181_01 


1997.02.03 


5000 


1.06 


1.02 


4.95 




d 9.1 + °- 2 7.0+°i 

— 0.2 1 —0.7 


P901&1 09 

IT iUlOJ \J£i 


1 QQ7 0.9 1 0. 


10500 


0.94 


0.91 


4.58 




do 9 +0.2 R q+0.7 
°- -0.2' "-^-o.r 


.r zUlol_Uo 


1 OQ7 no 1 7 


oZUU 


n on 
u.yu 


n c7 


4. Do 




do 7+O.2 f. ,+0.7 


P20056-01 


1997.04.05 


2200 


0.91 


0.87 


4.29 




F2UU5o_U2 


1997.04.10 


2100 


1.09 


1.03 


4.77 






P20056-03 


1997.04.11 


1800 


1.09 


1.02 


4.77 






P20056-05 


1997.04.15 


1800 


1.16 


1.06 


5.05 






1"! O /"V 1"! f / ' / \ / • 

P2UU5d_0d 


1997.04.17 


1300 


1.19 


1.11 


5.18 








1 OQ7 fi/l 1 

iyy / .U4.iy 


ziuu 


1.1 / 


l.U I 


A 7£ 
4. / 








1 QQ7 04 99 


2000 


1.10 


1.06 


4.42 






P20181_04 


1997.05.29 


9900 


0.60 


0.59 


3.02 






P20181.05 


1997.07.07 


10000 


0.24 


0.25 


1.48 






P20181.06 


1997.08.23 


11400 


0.72 


0.71 


1.84 






P20181.07 


1997.09.20 


10500 


0.95 


0.90 


4.49 






P20181.08 


1997.10.28 


10500 


0.62 


0.47 


3.21 






P40108.01 


1999.01.12 


12300 


0.63 


0.14 


0.20 




e < 0.2, < 0.2 


P40108.02 


1999.01.22 


13100 


0.53 


0.12 


0.20 




e < 0.2, < 0.2 


P40108.03 


1999.02.12 


6800 


0.47 


0.42 


1.84 


f 17.9, 17.2, 15.9 


e 4fi+ 01 6 4+ - 1 


P40108.04 


1999.03.03 


16600 


s 0.47 


0.45 


2.40 


9 17.6, 16.9, 15.6 


7+O.I fi 1+0.1 

°' '-0.1' D - 1 -o.i 


P40108.05 


1999.04.02 


9200 


90.49 


0.48 


2.75 




« 1+0-1 4 „+0.1 
°- ± -0.1' *-°-o.i 


P40108.06 


1999.04.22 


13100 


0.23 


0.23 


1.31 




e o , +0.0 9 9 +0.1 
°- 1 -0.0' z -^-o.i 


P40108.07 


1999.05.14 


9800 


0.08 


0.07 


0.38 




ei 4 +0.0 1 q+0.1 

1 -^-o.o> 1 -°-o.i 



"Ex posure times rounded to the closest 100 s; b PCA Fl ux: C H EXTE Flux, normalized to PCA; d 8.6GHz, 0.8 GHz data, from Wilms et 
al. ( 1999 ); e 8.6GHz, 4.8 GHz data, from Corbel et al. ( 200o| ); -^data from YALO telescope, 5-6 days after X-ray observation (/) and 
within ±lday of X-ray observations (<?), courtesy C. Bailyn and R. Jain; h GIS Flux 0.39 and 0.38 xlO -9 erg cm -2 s —1 , respectively 



ativistic beaming ( Shrader fc Titarchuk 1999|; jBcloborodov 



199£ 
2001 



Dove 



leavily ionized reflector mo d els flNayakshin 
f ayakshin fc Kallman 200l| ; pone 200 lj ); and mod- 



els in which the ADC and the accretion disc are physically 
separated. In the latter models, the ADC is typically repre- 
sented by a very hot, geometrically thick accretion disc in 
the central regions around the black hole, surrounded by a 
cold, geometrically thin and optically thick accretion disc. 



Model s ! of ftrte type are Lhe "riphere+diric 11 models of Dove 



et al. (1997b) whcie iL is assumed Lhal Lhe uiaLLei fluw- 



ing through the corona is efficiently converting its potential 
energy into radiation — and the Advection Dominated Ac- 
cretion Flows (ADAFs; e.g., Esin et al. 1997) — where most 
of the potential energy of the accreting matter is advected 
into the black hole. These latter models also postulate that a 
large fraction of the seed photons for Comptonization come 
from synchrotron radiation, due to the motion of the hot 
electrons in the magnetic field of the advective flow. 

A further complication of ADC models is the interface 
between the cold accretion disc and the ADC. In "slab-like" 
ADC models, the high ionization parameter at the ADC-disc 
interface leads to the formation of an ionized transition layer 
which can strongly aff e ct the reflectivity of the accretion disc 



flRoss fc Fabian 1993; |Nayakshin fc Dove 200 lj ; jNayakshin 
fc Kall jnan 200lj [Done 200l| ). The "sphere+disc" models, 
notably, fail to reproduce the observed iron fluorescence fea- 
tures, possibly due to the fact that any overlap between the 
accretion disc and the ADC is usually not modelled (Dove et 
al. 1998, although see Zdziarski et al. 1998). Finally, due to 
the closeness of the ADC to the central black hole, relativis- 
ts smearing of spectral features needs to be accounted for 



(Done & Zycki 1999), although many models have ignored 
such effects. 

Observationally, the situation is equally challenging. For 
"sphere+disc" models, radiatively efficient spherical mod- 
els have suggested coronal radii as small as ~ 30 GM/c 2 
(Wilms et al. 1999) . while ADAF m odels typically prefer 
radii j> 200 GM/c (Esin et al. 1997). Attempt s to measure 



the co r onal size mostly hav e been unsuccessful ( Nowak et al 



1999a 



Nowak et al. 1999c). On the other hand, attempts to 



model the observed broad band X-ray and 7-ray spectra 
with the ADC models have generally been quite successfu l 
(|Dove et al. 1998[ |Gierlihski et al. 1997|; |CierlMski et al. 199fi|; 



Zdziarski et al. 1998; Shrader fc Titarchuk 1999). Many fun- 



damentally different Comptonization models appear to give 
equally good fits to the data, thus making a decision between 
these models from a purely observational point very dim- 
cult. One potential discriminant among models has e merged 
with the suggestion of Zdziarski, Lubiriski, & Smith ( 199!: ) 
that there is an anti-correlation between the spectral hard- 
ness of hard state BHC spectra and the covering fraction of 
cold, reflecting material (the so-called T-fi/27r correlation', 
where F is the photon number flux power law spectral in- 
dex, i.e., F-y oc E~ r , and Q/2ty is the covering fraction of the 
cold reflector). Different power- law+reflection models which 
fit the same data equally well, have implied very different 
underlying fit parameters. (Contrast the fits of GX 339—4 
data presented by Wilms et al. 1999 to those presented by 
Revnivtsev, Gilfanov, & Churazov 2001.) Finally, we note 
that the "Comptonization parad igm" itself has recently been 



challenged, with Markoff et al. (2001) postulating that syn- 
chrotron radiation from the radio outflow/jet, now observed 
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in many BHC systems during the hard state, might well 
contribute a large fraction of the observed X-ray flux. 

1.2 Observations of GX 339 4 

The galactic black hole GX 339—4 is well-suited to perform 
a new study of the relative merits of different accretion disc 
models, as it is the only persistent sou rce that has been ob- 



Grebenev et al. 1991 




Miyamoto et al. 1991 




Mendez & van 


der Klis 1997 




Wilms et al. 1999|). Furthermore, the emis- 



sion from GX 339—4 is likely completely dominated by the 
accretion flow, all the way from radio through gamma-ray 
wavelengths, a nd it shows strong corre l ations among these 



energy bands (Hannikainen et al. 1998; Wilms et al. 1999 



Corbel et al. 2000; Fender 2001). In general, radio emission 



at 5i 10 mjy levels is present, is positively correlated with 
the X-rays in the hard state, and is quenched during the 
soft state ( 



Hannikainen et al. 1998 ; Fender et al. 1999 



Cor- 



bel et al. 2000] ). Brightness temperature arguments suggest 



that, during the hard sta te, the radio emittin g outflow must 
extend to > 10 12 5 cm ([Wilms et al. 1999|). For the hard 



state, the optical light from this system is likely dominated 
by the accretion flow onto the compact object as well, with 
a large fraction of the hard state optical flux possi bly being 
synchrotro n emission from the extended outflow ( [mamura 
et al. 1990; gteiman-Cameron et al. 1990| ). 



In this paper we apply a variety of spectral models 
and timing analyses to 22 observations of GX 339—4 ob- 
tained with the Rossi X-ray Timing Explorer (RXTE), 
some of which were also simultaneous with the Advanced 
Satellite for Cosmology and Astrophysics (ASCA), optical, 
and/or radio observations. The observation log, with fluxes 
in selected energy bands, is presented in Table |lj Assum- 
ing isotropic emission, a 4Mq compact object mass, and 
a distance of 4kpc, the 3-200 keV luminosities of the hard 
state observations listed in Table |l| correspond to 2 x 10 -3 - 
0.03 L Ed d. 

The observations were chosen for analysis as follows. 
Eight of these observatio ns (P20181) were originally ex- 



199E 



Nowak et al 



tensive ly analyzed by us (Wilms et al. 
1999b). We had found that those data, with the exception 
of the faintest observation, showed few differen ces among 
their variability properties (Nowak et al. 1999b), and that 



their spectra could be well-fit by either a sphere+disc coro- 
nal model, or by a p ower-law plus ionized reflection model 
(Wilms et al. 1999). For the latter spectral model, we 
did not find a hardness-reflection fraction anti-correlation, 
but instead we found a hardness-ionization parameter anti- 
correlation. This is in contrast to the results of Revnivtsev, 



Gilfanov, & Churazov (2001), who, analyzing the same data 
with a power law+reflection model (with gaussian smearing 
applied to mimic relativistic smearing), did find the sug- 
gested anti-correlation. We have chosen to reanalyze these 
observations, as well as seven others (P20056) that repre- 
sented the brightest and softest hard state obs ervat ions pre- 
sented by Revnivtsev, Gilfanov, & Churazov (2001). 

We present seven other observations (P40108) that are 
original to this work. All of these observations were per- 
formed simultaneously with radio observations (see Corbel 
et al. 2000), two were performed simultaneously with Ad- 
vanced Satellite for Astrophysics and Cosmology (ASCA) 



observations, and three were performed near-simultaneously 
with optical monitoring. An additional five observations ex- 
ist in this series (again, with a number of simultaneous radio, 
optical, and ASCA observations); however, these latter ob- 
servations, in contrast to the first seven, are too faint to 
perform useful timing analyses. These fainter observations 
will be discussed in a future work (Corbel et al. 2001; in 
preparation). 



1.3 Outline 

In §^ we discuss the spectral analysis in detail, concentrating 
on the results obtained with three different Comptonization 
models. Before discussing the analyses in detail, however, we 
begin with a comparison of the relative strengths and weak- 
nesses of each of these models. In §^ we describe the timing 
behaviour of GX 339—4, and correlate this behaviour with 
the results of our spectral analyses. Twenty out of twenty 
two of our observations occurred in the hard state. In §^| we 
describe o ur a nalysis of the two soft state observations. We 
assess in §5.2 the robustness of our fits by presenting ratios 
of the GX 339—4 data to observations of the Crab nebula 
and pulsar. We explore the r eali ty of the hardness-reflection 
fraction anti-correlation in 



5.3 



In §5.4, we present a new, 



simple coronal model that reproduces some of the coronal- 
temporal correlations that we find in §^. We discuss the 
implications of our results for ADAF models in §5.5, and we 
briefly discu ss alternative models to thermal Comptoniza- 
tion in § |5.6| We summarize our results in §^j. 

Our results rely on many observations with several 
satellites. In order not to deter the reader from the results 
of the analysis with the technicalities of the data extrac- 
tion, we present our data extraction strategies separately in 
Appendix [a|. 



2 SPECTRAL FITS 

2.1 The Blind Touching the Elephant: The 

Relative Merits of Various Coronal Models 

Here we d iscuss the Co mpton ization codes of Poutanen & 
Svens son (|l996| ), Coppi ( |l999| ), and Dove, Wilms & Begel- 
man (1997a). Although each has been used separately to fit 



various BHC spectra, this work is the first to attempt to 
systematically compare several of them. All three will be 
used in fits of the GX 339—4 data. Each model has unique 
attributes and considers one or more facets of the Comp- 
tonization problem more 'correctly' or 'completely' than the 
others; however, none are truly fully se lf-con sistent models. 

The code of Poutanen & Svensson (1996) (compps in the 
X-ray spectral fitting package, XSPEC; Arnaud 1996) is es- 
sentially a 'one-zone' model in that it considers injection of 
(either blackbody or multi-temperature disc blackbody; Mit- 
suda 1994) seed photons with a prescribed flux and a fixed 
geometry into a corona of uniform optical depth and temper- 
ature and/or uniform distribution of non-thermal electrons. 
The corona also has a fixed geometry, although a variety of 
geometries are considered. A fraction of the escaping pho- 
tons can be reflected off of cold, partially ionized material; 
however, this fraction is not self-consistently determined, nor 
are any of the reflected photons further reprocessed by the 
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corona. Smearing of the reflection due to relativistic disc 
motion and gravitational effects, however, can be applied. 
The Fe edge is modelled in the reflection, but the Fe line is 
not. Furthermore, the temperature of the corona and its 
(photon-photon collision produced) electron-positron pair 
optical depth are not self-consistently calculated, and must 
be checked a posteriori. 



The model of Coppi (1999), available from that author 



as the XSPEC model eqpair, does calculate a self-consistent 
coronal temperature and pair optical depth by using the 
coronal compactness (proportional to energy released in the 
corona divided by its radius) as the fundamental fit param- 
eter. Furthermore, it allows one to specify both a thermal 
and non-thermal compactness for the corona, as well as a 
compactness for the (blackbody) seed photons. A spheri- 
cal corona with seed photons distributed according to the 
diffusion equation is modelled; therefore, all the seed pho- 
tons pass through the corona. Reflection is implemented in a 
similar manner to the compps model, but smearing of the re- 
flection is not considered. We will use this model extensively 
in the following sections, and we will present coronal com- 
pactness, £ c , instead of coronal temperature, kT a , coronal 
2/-parameter, or power-law photon index, V, as a fundamen- 
tal fit parameter. 

The coronal model of Dove, Wilms, & Begelman 



(1997a), available from those authors as the XSPEC model 
kotelp, utilizes Monte Carlo spectra calculations that are 
then stored in an interpolation grid. For a given seed pho- 
ton spectrum and geometry and a given coronal geometry 
and initial optical depth, the non-uniform coronal tempera- 
ture distribution, pair optical depth, energy balance between 
corona and disc, and reflection fraction — with subsequent 
re-Comptonizations/re-reflections of this component — are 
all self-consistently calculated. In this work we sh all us e the 



'sphere+disc' geometry described by Dove et al. (1996), i.e. 
a multi-temperature disc blackbody seed photon distribu- 
tion with kT(R) ex R~ 3 ^ and maximum temperature of 
kT = 150 eV emanating from an outer thin disc surround- 
ing a central spherical corona. (Other geometries and seed 
photon distributions are possible.) Again, the fundamental 
parameter is the coronal compactness, £ c , relative to the disc 
compactness (fixed to 1), not the average coronal temper- 
ature. For £ c Si 1, the intrinsic energy generation of the 
disc dominates the system flux. For £ c J; 1, the seed pho- 
ton spect rum is dominated by hard X-rays reprocessed by 
the disc (Dove et al. 1997b). The Fe line and edge are both 
calculated as part of the reflection component; however, in 
current implementations no relativistic smearing and only 
solar abundances are considered. 

As is evident from the above capsule summaries, the 
various models do represent to some extent the "blind man 
touching the elephant", each feeling out a different portion 
of the puzzle. Even though the compps model represents a 
'one zone' corona, it is still the most flexible in terms of ge- 
ometry of both the corona and the properties of the reflector. 
(Given sufficient computer resources to run grids of models 
kotelp can be equally flexible.) The compps model, however, 
lacks the self-consistent energy balance provided by eqpair 
or kotelp. The eqpair model has self-consistent energy and 
pair balance, but lacks the flexible geometry of the other two 
models. The kotelp model self-consistently determines the 
coronal energy balance, the seed photon flux (i.e., intrinsic 



disc flux vs. reprocessed flux), and the reflected component, 
but is the least flexible of the models in terms of currently 
available fit parameters. 

2.2 Considerations for Model Fitting 

A complication in utilizing the above models is that de- 
generacies are common, especially as regards the issue of 
reflection. Models of the Fe line/edge region are strongly 
affected by assumptions concerning the Fe abundance, the 
shape of the seed photon and incident hard spectrum, the 
reflection fraction, the reflector's ionization parameter, the 
coronal geometry and outflow properties (e.g., Beloborodov 
1999), and the reflector inclination. All but the first of this 
list could easily vary from observation to observation. (The 
reflector inclination could vary due to disc warping effects; 
Pringle 1996, Maloney & Begelman 1997.) Indeed, we have 
found adequate fits for the same observation with reflec- 
tion fractions ranging from to 1, Fe line equivalent widths 
ranging from 50-400 eV, coronal temperatures ranging from 
20-200 keV and coronal optical depths ranging from 0.1-3. 

Another consideration perhaps unfamiliar to those who 
do not practice X-ray astronomy and perhaps too readily ac- 
cepted by those who do, is that model fits often have reduced 
\ 2 substantially below unity. This is in part due to the adop- 
tion of systematic errors (see the Appendix), and indicates 
that we are fitting some systematic features of the satellite 
response matrices. Fits to the (presumed power-law) Crab 
nebula plus pulsar clearly indicate the presence of systematic 
errors in the PCA response. Furthermore, without system- 
atic errors, the low energy PCA data would dominate the 
fits, while the high energy HEXTE data, useful for mea- 
suring coronal temperatures and optical depths, would add 
very little statistical weight. 

Obtaining a reduced \ 2 substantially below one indi- 
cates that an individual observation is over-parameterized. 
In the fits described below, we try to minimize any over- 
parameterization by freezing numerous parameters, espe- 
cially in the line region, and we have tended to adopt fairly 
small systematic errors in high signal-to-noise regions of the 
spectrum (e.g., below 4keV). The reduced \ 2 f° r the fits we 
present are always > 0.6 (cf. Revnivtsev, Gilfanov, & Chu- 
razov 2001). We would consider this too small if we were 
fitting a single observation; however, we are concerned with 
the relative trends of the fit parameters. Here we consider 
20 separate, but uniformly measured and uniformly fit hard 
state observations. For example, small differences in the line 
region, even when some fraction of this region represents 
systematic efects, can be very statistically significant given 
the large effective area of the PCA. Differences in fit param- 
eters, therefore, can be meaningful. Even then we endeavor 
not to over-interpret the parameter trends (see, for example, 
the discussion of Fig. ^| ). To give the reader some indi ca- 



Al 



tion of the role played by systematic errors, in Table 
present a set of fits for observation P20181_01 (italicized), 
where we adopt no systematic errors. The error bars for pa- 
rameters associated with narrow features (e.g., line widths) 
decrease somewhat. Overall, however, the "bestfit" param- 
eters are not substantially changed, although the \ 2 has 
greatly increased. 

Given the above concerns about model degeneracies and 
systematic errors, we explicitly examine the ratios of our 
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10 20 50 

Energy (keV) 

Figure 1. Unfolded energy spectrum of GX 339-4 (P40108.03) 
with the eqpair Comptonization model fitted to the PCA data 
only. The HEXTE data, assuming the model fit to the PCA data, 
is also shown assuming a relative normalization constant of 0.9 
(see Appendix) for both the HEXTE A and B clusters as com- 
pared to the PCA. (Here error bars arc l-cr.) 



data to those from the Crab nebula and pulsar. We discuss 
the implicatio ns o f these comparisons for the interpretation 
of our fits in 



5.2 



We further discuss the relative parameter 



trends in the sections below. 



2.3 PCA Fits: eqpair+gauss 



As discussed by Wilms et al. (199!:), there typically has been 
a difference between the best-fit photon indices for PCA 
and HEXTE observations of the Crab pulsar and nebula 
of AF J; 0.1. Recent PCA response matrices have slightly 
decreased this discrepancy; however, differences still remain. 
Partly for this reason, numerous authors have chosen not to 
include the HEXTE data in their fits of GBHC spectra (e.g., 
Revnivtsev, Gilfanov, & Churazov 2001). Here we briefly 
explore whether this is justified. 

For our fits to PCA data, we have applied the eqpair 
Comptonization model with an additional gaussian line pro- 
file with a peak energy fixed at 6.4 keV. We further have 
fixed the blackbody seed photon compactness to be 1, the 
reflection to be from cold, neutral material with solar abun- 
dances, the angle of the reflecting medium with respect to 
our line of sight to be 45° , and the compactness of the non- 
thermal electrons to be 0. The fit parameters are therefore 
the compactness of the thermal electrons, the coronal optical 
depth, the seed photon blackbody temperature, the over- 
all normalization, the line width, and the line amplitude. 
Folding in the systematic uncertainties described in the Ap- 
pendix, these models yield extremely good fits to all the 
hard state observations, with typical reduced \ 2 ~ 0.3. The 
fitted lines have widths of a « 1 keV and equivalent widths 
ranging from 100-300 eV. The reflection fractions range from 
0/27T ~ 0.3-0.6. The seed photon temperature is typically 
ps 600 eV, the coronal optical depth is typically » 2, and 
the coronal compactness is typically £ c ~ 1-7, which here 
implies coronal temperatures of tens of keV. 

Fig. [j], which presents the unfolded spectrum for obser- 
vation P40108_03, shows that these fits fail utterly to re- 



produce the HEXTE data, even accounting for likely sys- 
tematic differences between PCA and HEXTE. The 'un- 
folded spectrum' shows the best fitting model, here plotted 
as E 2 x Fj(E), multiplied by the ratio of the data in detector 
space to the spectral model folded through the experimental 
response function of the detectors. It is important to note in 
Fig. [l] that the apparent visual significance of narrow spec- 
tral features should not be taken too strongly. An 'unfolded 
spectrum' has a tendency to exaggerate the assumed spec- 
tral features, especially those with widths comparable to or 
narrower than the instrumental resolution. In this instance, 
however, the unfolded spectrum does highlight the complete 
inadequacy of the extremely well-fitting PCA model in de- 
scribing the broad-band HEXTE data. 



2.4 PCA+HEXTE Fits: kotelp Models 

We fit the simultaneous PCA and HEXTE data with the 
kotelp mod el wit h a 'sphere+disc' geometry as described by 
Dove et al. ( 199£ ) . Here we choose to average the spectrum 

These models do 



over inclination angles /i = cos 9 = 0.4-0.1 
not adequately fit the data, and typically produce reduced 
X 2 > 3, primarily due to large, broad residuals in the Fe 
line/edge region, as one might expect for a relativistically 
broadened line. (We previously had noted this for kotelp 
fits to RXTE data of Cyg X-l; Dove et al. 1998.) 

In order to mimic such a line for the kotelp model, we 
performed a series of fits where we added a narrow, fixed en- 
ergy and negative amplitude gaussian function to remove the 
« 30 eV equivalent width (EW) line produced by the code. 
We then added a broad gaussian function with peak energy 
fixed at 6.4 keV, and refit the models. Fitted parameters 
are presented in Table 



Al 



in the Appendix. These models 
produced markedly improved fits, with reduced \ 2 ranging 
from 0.7-2.5, with an average reduced \ 2 °f 1-6- The three 
faintest observations fit best. Most observations show neg- 
ative residuals, as might be expected from an unmodelled 
edge, near lOkeV. A large fraction of x 2 > however, comes 
from the 3-4 keV region, where 1999 observations show pos- 
itive residuals, i.e., a soft excess, and 1997 observations show 
negative residuals. 

We note that if we instead fit near edge-on kotelp mod- 
els, the P20056 fits improve by A\ 2 ~ 10, since the contri- 
bution of the seed photon spectrum is reduced compared to 
the Comptonized spectrum, and hence the negative 3-4 keV 
residuals are reduced. The P20181 observations fit approx- 
imately the same, but the 1999 fits become much worse as 
the 3-4 keV excess increases. If we instead choose a nearly 
face-on kotelp model, as the seed photon spectrum is now 
emphasized compared to the Comptonized spectrum, the 
positive 3-4 keV residuals are then reduced. Such fits to the 
1999 observations improve by Ax 2 » 30. The fits to the 1997 
observations, however, beco me much worse. We discuss these 
points further in §2.5, §2.6, and § [5.2| . 

Fitted compactnesses range from 1-7, and fitted op- 
tical depths range from 1—1.6. These parameters translate 
to average coronal temperatures of 90-170 keV. (Note, how- 
ever, that face-on kotelp models produce lower compact- 
nesses and higher optica l de pths, resulting in a corona with a 
lower temperature; see §2.6 ) Equivalent widths of the added 
broad line range from 110-220 eV, and widths of the lines 
range from a ~ 0.4-0.9. 
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Figure 2. Results of the kotelp model fit to joint PCA and HEXTE data of GX 339—4. Top: Coronal compactness vs. PCA flux in the 
3-9 kcV band. Bottom: Equivalent width of the Fe line vs. coronal compactness. Throughout this work, unless otherwise noted, diamonds 
represent the 1997 hard state observations previously discussed by Wilms et al. (1999). Triangles represent 1997 hard state observations 
discussed by Revnivtsev, Gilfanov, &: Churazov (2001). Squares represent 1999 hard state observations. Circles represent the first 1999 
hard state observed after GX 339—4 returned from a protracted soft state. The dotted vertical line represents the lowest 3— 9keV PCA 
flux observed for the soft state observations discussed here. The dotted error bars are the 90% confidence limit (i.e., Ax 2 = 2.71 for one 
interesting parameter), and the solid error bars are equal to the 90% confidence level error bars divided by V2.71. 



Two trends stand out from these fits. As Fig. |^ shows, 
there is a strong anti-correlation between compactness and 
3-9 keV PCA flux for the 1997 hard state observations. In 
1999, GX 339—4 is seen to have returned from the soft state 
with an initially very soft spectrum that hardened as the 
source faded, before beginning to soften slightly at very low 
flux levels. As previously noted with Gmga data (Miyamoto 
et al. 1995), there is clear evidence for hysteretic behaviour. 
In 1997, GX 339-4 achieved a higher 3-9 keV flux, with- 
out entering the soft state, than the lowest 3-9 keV soft 
state flux measured shortly prior to the return to the hard 
state in 1999. The 'overlap region' between the hard and 
soft states appears to be in the 3-9 keV PCA flux range 
of ~ 0.5-1.2 x 10~ 9 ergs cm 2 s _1 . Of course, the issue of 
the ordering of the bolometric luminosities is less clear, es- 



pecially since we have not measured the flux below 3keV, 
which likely accounts for the largest fraction of the soft state 
flux. The observations represented by triangles in these fig- 
ures are among the highest flux levels ever detected for the 
hard state of GX 339—4, whereas other 'high state' and 'very 
high state' obse rvations have exhibited substantially great er 
3-9 keV fluxes ( feelloni et al. 1999[ |Miyamoto et al. 199 1| ). 

The second correlation that stands out, also shown in 
Fig. ^, is the anti-correlation between the equivalent width of 
the broad line and the coronal compactness. In previous fits 
to a fraction of the 1997 hard state observations (diamonds; 
Wilms et al. 1999), utilizing PCA data only, such a corre- 
lation was not detected. Aside from not including HEXTE 
data, those previous fits utilized an earlier version of the 
PCA response matrix with larger systematic uncertainties. 
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Here to some extent the Fe line equivalent width might be 
serving as a reasonable proxy for reflection fraction. If this 
is the case, however, then whereas a correlation is seen prior 
to the 1997 transition to the soft state, the correlation is 
apparently absent after the 1999 return to the hard state. 

We considered two other combinations of models with 
the basic kotelp model in order to assess the robustness of 
the above correlations. Here we shall discuss only the high- 
lights of the results. For both of these models, the line pro- 
duced by the kotelp code was again subtracted as above. I n 
the first case, we added the XSPEC laor model ( [Laor 199l| ), 
which gives the line profile for a relativistically broadened 
line from a disc rotating about a Kerr black holeQ. These 
models improved the reduced x 2 of the fits to 0.7-2.1, with 
an average reduced % 2 of 1.2. (The observation immediately 
following the return to the hard state, P40108_03, produced 
the worst fit, again due to a soft excess at ^ 4keV.) The 
line energy, which was allowed to be a free parameter, ranged 
from 6.5-6.8 keV, and the line equivalent widths ranged from 
150-350 eV. There essentially was no change from any of the 
trends and correlations discussed above. 

For the second model, in addition to using the laor 
model to describe the line region, we also added an ioniza- 
tion edge at ~ 8keV. This further improved the fits to a 
reduced \ 2 ranging from 0.7-1.4, with an average of 0.9. 
(As before, observation P40108_03 was the most poorly fit 
due to a soft excess at energies Js 4keV.) The line energies 
were as above and their equivalent widths ranged from 100- 
270 eV. The additional ionization edge had an energy that 
ranged from 8.0-9.5 keV and that also was positively cor- 
related with the laor line energy. The edge optical depth 
ranged from r = 0.02-0.07. Although the line equivalent 
width correlations discussed above were slightly weakened, 
the trends shown in Fig. ^ again remain basically unaltered. 
Taking the above results as a potential signature of reflection 
changes, we turned to fits with the eqpair model. 



2.5 PCA+HEXTE Fits: eqpair Models 

For the joint fits to the PCA and HEXTE data d iscu ssed 
here, we used the exact same eqpair model as in §2.3, ex- 



cept that we instead fixed the Fe abundance to be 4 times 
solar, which is more consistent with the large line equiv- 
alent widths but low reflection fractions. (Lower or higher 
abundances, however, produced only slightly worse fits.) We 
found acceptable fits to the joi nt da ta with parameters very 
similar to the PCA only fits ({2.3); however, a completely 
different set of parameters described the joint data slightly 
better. These models, with a reduced x 2 « 0.6-1.0, had seed 
photon temperatures in the range of 30-100 eV, coronal com- 
pactnesses in the range of 5-14 (yielding coronal tempera- 
tures of ~ 200 keV, given the typically fitted optical depth of 
0.1-1), and reflection fractions of 0.1-0.5. Line widths were 
broad, and line equivalent widths were ~ 80-240 eV. Since 
these are considered to be the more "usual parameters" (see 



* To some extent, the broad red tail of the laor line can off- 
set differences between the HEXTE and systematically softer 
PCA detectors. A hard power law extrapolated backward from 
the HEXTE energy band to the PCA energy band could leave a 
soft excess, even if in reality no such excess existed. 



Poutanen, Krolik, & Ryde 1997) , we present these results in 
Figs. and give the parameters in Table Al. 

As shown in Fig. |§| the same basic pattern of com- 
pactness vs. 3-9 keV flux was seen, i.e., higher fluxes yielded 
softer spectra, and GX 339—4 returned in 1999 from the soft 
state with an initially soft spectrum. The variation, how- 
ever, is less pronounced, especially on the return from the 
soft state in 1999. This is because for the eqpair model 
the seed photon temperature is not a fixed quantity. The 
relatively soft spectra for the first three 1999 hard state 
observations (P40108_03-_05) are modelled with a combi- 
nation of both low compactness and low seed photon tem- 
perature (« 30 eV). We note also that these three spectra 
indicate coronae that are pair-dominated, as the fitted ini- 
tial electron optical depth was r cs = 0.01 (we set this value 
as a lower limit to the fits), but the total optical depths 
were r ~ 0.3. The error bars on the initial electron op- 
tical depths, however, were large, and non-pair-dominated 
plasmas are allowed. The soft positive residuals seen for ob- 
servation P40108_03 with the kotelp model fits were absent 
with the eqpair fit. 

Reflection fraction and Fe line correlations are readily 
noticed in Figs. ^| and ^. The most striking pattern is the re- 
flection fraction, fl/2n, vs. the 3-9 keV PCA flux. There is a 
positive correlation of fitted reflection fraction with flux, ex- 
cept in the overlap/hysteretic region between the hard and 
soft states. There the reflection fraction is quasi-independent 
of flux, except at the very highest flux levels (triangles). 
Correlating the reflection fraction with compactness is less 
clear. Ignoring the P20056 observations (triangles), and the 
3 lowest flux observations (P20181_05, P40108_06-_07), no 
correlation would be seen between reflection fraction and 



compactness, as discussed by Wilms et al. ( 1999 ). Only by 
including the highest flux (triangles) and lowest flux obser- 
vations does a possible anti-correlation between compact- 
ness and reflection fraction appear. Furthermore, the faintest 
three observations, taken by themselves, indicate exactly the 
opposite trend. 

The Fe line equivalent width, as for the reflection 
fraction, is anti-correlated with compactness if one con- 
siders all the observations, but does not show this anti- 
correlation if one ignores the brightest (triangles) and 
faintest (P20181_05) observation (see Wilms et al. 1999). 
This is also seen by correlating the Fe line equivalent width 
with the reflection fraction. Prior to the 1997 transition to 
the soft state, reflection fraction and Fe line equivalent width 
are correlated, after the 1999 return to the hard state, there 
is a lower limit to the Fe line equivalent width of ~ 150 eV. 
This is com para ble to what we found for the kotelp fits 
discussed in §2.4, except that there the 1999 post-transition 
equivalent widths were consistent with the lowest equivalent 
width measured prior to the 1997 soft state transition. 

There also is an apparent anti-correlation between phys- 
ical line width, a, and coronal compactness, as shown in 
Fig. P This is similar to the anti-correlation between spec- 
tral hardness and 'smearin g wid th' claimed by Revnivt- 
sev, Gilfanov, & Churazov (2001). The low reduced \ 2 of 
our fits, however, indicate that we may already be over- 
parameterizing the data. Furthermore, this anti-correlation 
is greatly weakened if we fit the brightest observations (tri- 
angles) with fixed line energies of 6.7 keV instead. The line 
widths then decrease by Act m 0.2 keV. 
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Figure 3. Results of the eqpair model fit to joint PCA and HEXTE data of GX 339—4. Left: Coronal compactness vs. PCA flux in the 
3-9 keV band. Middle: Reflection fraction vs. PCA flux in the 3-9 keV band. Right: Reflection fraction vs. coronal compactness. 
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Figure 4. Results of the eqpair model fit to joint PCA and HEXTE data of GX 339-4. Left: Equivalent width of the Fe line vs. 
reflection fraction, fi/27r. Dashed lines show: the expected trend if the line equivalent width is proportional to reflection fraction (with 
intercept), and the mean line equivalent width for the P40108 observations. Middle: Equivalent width of the Fe line vs. PCA flux in the 
3-9 keV band. Right: Equivalent width of the Fe line vs. coronal compactness. 
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Figure 5. Width of the Fe line vs. the coronal compactness, for the eqpair model fit to joint PCA and HEXTE data of GX 339—4. 
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2.6 ASCA Observations 

Two of the RXTE observations discussed in this work, 
P40108JD4 and P40108_05, were obtained simultaneously 
with ASCA. Compared to RXTE, ASCA had a nar- 
rower band pass (effectively 1-10 keV for observations 
obtained as late as 1999), an approximately 4 times 
smaller effective area, but a much higher spectral resolu- 
tion 0.2 keV for these observations) and a lower back- 
ground due to the nature of its instruments and its ar- 
cminute spatial resolution. Furthermore, as it is an indepen- 
dently calibrated instrument (see the mission web pages at 



tittp : / /heasarc . gsf c .nasa.gov/docs/ asca/ ascagof . html 

for a discussion of ASCA calibration issues), it provides a 
useful check of the RXTE results. 

Analyzed independently of RXTE, these two ASCA ob- 
servations can be fitted wit h mod els very similar to those 

employed by Wilms et al. |l \ to study three archival 

ASCA observations of GX 339—4. Specifically, these obser- 
vations are well fit by a multi-temperature disc blackbody 
spectrum with peak temperature of ~ 0.15 keV, a broken 
power law with a photon index ~ 1.8 below 4keV and ~ 1.6 
above, and a gaussian line at 6.4keV[J] These observations 
prefer Fe lines with equivalent widths of ~ 220 eV and that 
are broad (a » 0.8 keV), although the exact parameter val- 
ues for the line do depend upon the fitted continuum model. 

One might hope that including the RXTE data would 
remove enough of the ambiguities in the continuum model so 
as to further constrain the line model. In Fig. |^ we plot the 
residuals (detector space data divided by the model folded 
through the response matrices of the detectors) for an ab- 
sorbed power law (jV H fixed to 6 x 10 21 cm~ 2 ) fit to 1-22 keV 
data (1-10 keV in CIS, 2-10 keV in SIS, and 3-22 keV in 
PC A). Although both ASCA and PC A agree on the pres- 
ence of a strong, broad line-like feature near 6.4 keV, there 
are clear systematic differences among the CIS, SIS, PC A, 
and HEXTE detectors. As we note in the Appendix, PCA 
and HEXTE yield different power law indices for observa- 
tions of the Crab nebula and pulsar, with HEXTE giving a 
systematically harder index. The ASCA detectors were cali- 
brated to yield Crab indices comparably hard to that found 
by HEXTE. Furthermore, both the SIS and CIS detectors 
show a break near 4keV that is absent in the PCA data, as 
well as additional spectral structure near 1 and 2keV. 

Ignoring the PCA data, data from the CIS, SIS, and 
HEXTE detectors can be reasonably well-fit with the same 
disc blackbody, broken power law, plus gaussian line models 
mentioned above, but with an additional caveat. The power 
law above 4keV (uniform from ASCA through HEXTE en- 
ergy bands) is exponentially cutoff in the HEXTE bands at 
energies only i; 80 keV (i.e., the highecut model in XSPEC; 
see also Wilms et al. 1999). The reduced \ 2 f° r such models 
are m 1.3, with the greatest deviations coming from the 8- 



t Note that, consistent with the analy sis of A SCA observati ons o f 
Cyg X-l presented by Ebisawa et al. (1997), Wilms et al. (L999) 
searched for fits with narrow (a ~ 0.1 keV) Fe lines. The fitted 
lines had equivalent widths m 40 eV, and the fits yielded reduced 
Y 2 ~ 1 for the faintest observations, but 1.4 for the brightest. As 
compared to the ASCA observati ons di scussed here, the observa- 
tions discussed by Wilms et al. ( 1999[ ) had 3— 9keV fluxes that 
were approximately 0.3, 0.5, and 1.8 times as great. 



10 keV and 17-25 keV regions, where the model over-predicts 
the flux. An even better fit to just the ASCA+HEXTE data, 
however, can be obtained with a near face-on kotelp model 
with the addition of a dust scattering halo, as is kn own t o 
exist in this system from the work of Predhel et al. ( 1991 ) , 
and references therein. We use the dust model from XSPEC, 
except here we add a component by fitting in XSPEC: ab- 
sorption x (model + [model — dust x model]). This accounts 
for scattering of X-rays into our line of sight. The extra halo 
component has two parameters: /halo, which is essentially 
the amplitude of the scattered component, and S^aio, which 
is related to size of the halo and effectively determines the 
energy ranges scattered by the halo. Residuals for these fits 
(allowing the JVh column to be a free parameter) are pre- 
sented in Fig. [?|, and the parameters are listed in Table 



A2 



These fits to the ASCA+HEXTE data are extremely 
good, with reduced \ 2 ~ 1 for P40108.04. (The reduced \ 2 
for P40108J35 is larger predominantly due to some disagree- 
ment between the SIS and CIS detectors in the 8-10 keV 
range.) The fitted Ah column (chosen to be the same for 
both the point source and dust spectrum) is in good agree- 
ment with prior measurements as well as considerations from 
optical extinction and 21 cm measurements (see the discus- 
sion and references of Zdziarski et al. 1998). The fractional 
contribution of the dust halo is ~ 10% at 1 keV, in agree- 
ment with the results of Predehl et al. (1991). The break in 



the ASCA spectra near 4keV arises from a combination of 
the red tail of the Fe line, the curvature of the Comptonized 
spectrum towards low energy, the prominent blue tail of the 
seed photons due to the face-on disc, and the dust halo. 
We note that these face-on kotelp models, as compared to 
the models of §2.4, have slightly increased coronal compact- 
nesses, but greatly increased coronal optical depths. This 
yields average coronal temperatures of ~ 50 keV, as opposed 
to the : 



2.4 



150 keV for the models of 1 
If we add the PCA data to the fits, the neutral hydrogen 
column and dust parameters do not greatly change, but the 
coronal compactness increases and the optical depth drops. 
The overall fit in terms of reduced \ 2 i s n °t as good as 
for the ASCA+HEXTE fits alone, as is clear from Fig. 0. 
The detectors, as for Fig. [], again show clear systematic 
differences. However, the level of residuals shown in Fig. ^, 
i.e., < 1-5% for the CIS, SIS, and PCA detectors, are as 
low one can reasonably expect to achieve given the differing 
systematic uncertainties of these detectors. 



3 TIMING ANALYSIS 



As discussed in { 1.2, these observations were selected partly 
because with them we were able to perform timing as well 
as spectral analyses. Timing analyses of the P20181 obser- 
vations p reviou sly have been prese nted by Nowak, Wilms, 
& Dove ( 1999b ) and Nowak ( 2000 ). Fourier power spectral 
densities (PSD) only of the P20181 and P20056 observations 
previously have been presented by Revnivtsev, Gilfanov, & 
Churazov (2001). The PSD fits p resent ed below most closely 
follow the discussion by Nowak (2000). 

Following the work of Wijnands fe van der Klis ( 199!: ) 
and Psaltis, Belloni, & van der Klis jl999|), Nowak (200C) 



suggested that the best representation of the hard state PSD 
of GBHC in general, and GX 339—4 in specific, was the 
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Figure 6. Residuals for an absorbed power law model fit simultaneously to the SIS, GIS, and PCA data (P40108-04). Circles and 
crosses are the GIS and SIS data, respectively. Diamonds are the PCA data. Stars are the HEXTE data (not included in the fit), with 
normalization constants chosen to match the PCA residuals. (Error bars here are 1-cr errors.) 
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Figure 7. Residuals for corona plus dust scattering halo models (see text) fit to: PCA and HEXTE data (top); GIS, SIS, and HEXTE 
data (middle); and GIS, SIS, PCA, and HEXTE data (bottom). (Error bars here are 1-cr errors.) 



s um o f only a few, broad, quasi-periodic features. Nowak 
(2000) showed that a composite PSD of the seven bright- 
est P20181 observations (which exhibit relatively uniform 
variability properties; see Nowak, Wilms, & Dove 1999, and 
the discussion below), could be described by a zero fre- 
quency centered Lorentzian function with a cutoff frequency 
at 0.1 Hz, plus additional quasi-periodic oscillations (QPO) 
at frequencies 0.34 Hz (and its harmonic), 2.5 Hz, and 18 Hz. 
Follow ing the suggestion of Psaltis, Belloni, & van der Klis 
(1999), the latter three frequencies might be the analogues 
of the 'horizontal branch', 'lower kiloHertz', and 'upper kilo- 
Hertz' QPO features seen in neutron star Z-sources. (For a 
description of the Z-source features, see van der Klis 1996.) 



Here we considered each observation individually; there- 
fore, we were unable to fit any features with peak frequencies 
<; 10 Hz. Furthermore, we chose to fit a zero frequency cen- 
tered Lorentzian function, 



P(f) = 



(1) 



- , (///c) 2 

to the lowest 0.1 Hz) and highest (J; 1Hz) frequency 
features. For the middle frequencies (/ « 0.3 Hz), we fit a 
QPO function of the form 



R 2 Qfo 



(2) 



/o 2 + W-/o) 2 ' 
where /o is the resonant frequency and Q is the quality 
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Figure 8. Fits to several PSD of GX 339—4. Solid line shows the best-fit model, while the dashed lines show the individual components of 
the model. Left: Observation P40108.05 in the 0-5 keV band. Middle: Observation P40108.05 in the 10-20 keV band. Right: Observation 
P20056-05 in the 10— 20keV band. (Error bars here are 1-cr errors.) 



factor, with Q fn A///o being the fractional width of the 
QPO. (For Q C 1, the maximum of P(f) occurs at / 3> 
/o. Partly for this reason, we instead chose to fit a zero 
frequency centered Lorentzian function at high frequencies, 
where the signal-to-noise was low. A number of observations, 
most notably P20056_05, were better fit, however, with a 
high frequency QPO feature. We fit multiple QPO features 
in the 0.3 Hz range, so long as the x 2 of the fits improved 
by more than 15 for each additional feature. Most PSDs 
required two such features, although several required three 
(Fig. §). 

Observation P40108_03 required only two features: the 
zero frequency centered Lorentzian function and a second 
feature which we identified as a QPO. As for all the other 
PSDs (with the exception of P40108.07; see below), this 
'QPO' represented a distinct maximum in the plot of fre- 
quency times the PSD power (e.g., Fig. ||). P40108_07 also 
showed only two features, but they both were broad, they 
were widely separated in frequency, and they were of com- 
parable amplitude in plots of frequency times PSD power. 
We therefore identified both of these features with the zero 
frequency centered Lorentzian functions. 

Fig. |§] shows several examples of the PSD fits, and pa- 
rameters fo r th e high energy band fits are presented in Ta- 
bles A4 and A5. In general, the root mean square (rms) vari- 



ability was on the order of 30%, with the soft energy band 
showing slightly higher rms variability than the high energy 
band. Decreased variability at high en ergy has been argued 
by Revnivtsev, Gilfanov, & Churazov (2001), via "frequency 
resolved spectroscopy" , as indicative of the "reflection spec- 
trum" having little variability. However, the PSD retains no 
phase information (which, in reality, could be a complicated 
mixture of components; see Nowak 2000). If a highly vari- 
able "power law" is 180° out of phase with a highly variable 
"reflection component" (as suggested for the u F-Q/2n" cor- 
relation; Zdziarski, Lubinski, & Smith 1999) , low variability 
would result at high energy. 

There is some trend for the PSD fit components to show 
lower rms variability with higher frequency. This trend is 
most pronounced for the highest frequency fit component, as 
we show in Fig. [j. We also show in Fig. ^the fitted frequen- 
cies plotted against the measured 3-9 keV flux. Whereas 



there is a large degree of variation among the low flux fre- 
quencies, there is relatively less frequency variation in the 
'hysteretic region' above the lowest measured soft state 3- 
9 keV flux. The P20056 observations at very high flux show 
a slight upward trend in frequency. Observations P20181_05, 
P40108_06, and P40108_07, all at fairly low fluxes, all exhibit 
frequencies markedly lower than those exhibited by observa- 
tions in the 'hysteretic region'. Observation P40108_03, the 
first observation after the return to the hard state in 1999, 
although exhibiting a relatively lower variability amplitude, 
shows a sharp spike upward in characteristic PSD frequen- 



Note that although there are clear systematic trends, 
there is not a very large variation in the characteristic fre- 
quency of a given PSD feature. For example, the cutoff fre- 
quencies of the Lorentzian functions fit at high frequencies 
vary b y less than a factor of 3 overall, di Matteo & Psaltis 
( [L999D have argued that the limited range of frequency varia- 
tions exhibited by GBHC argue for a comparably small (less 
than a factor of 2) variation in the (cylindrical) radial size 
scale of any central corona or in the size scale of any "transi- 
tion radius" between an advection and non-advection domi- 
nated flow. Assuming that the highest observed frequency is 
generated in the outer, non-advection dominated disc (e.g., 
Psaltis & Norman 2000), 1.8 Hz implies a transition radius 
270 GM/c 2 , for a 4Mq central object. However, a fea- 
ture wit h a fr equency as high as 20 Hz, such as discussed by 
Nowak ( 200C| ) (which would not be detectable within these 
individual observations) would imply a radius ^ 60 GM/c 2 . 
The rise of the characteristic PSD frequencies with increas- 
ing flux, followed by a plateau at high fluxes, is very rem- 
iniscent of the behaviour seen in the neutron star source 



4U 1820-30 flZhang et al. 1998; ; [Kaaret et al. 199J ). Such a 
plateau has been argued as evidence for the accretion disc 
extending inward all the way to its m arginally stable orbit 



flZhang et al. 199$ |Kaaret et al. 199§ ) 



In Fig. ^ we also show these frequencies plotted vs. 
the coronal compactness from the kotelp+gauss fits of Ta- 
There is a clear trend for harder spectra to show 



ble Al 



systematically lower frequencies for all of the characteris- 
tic PSD feat ures. Such a trend was noted by di Ma tteo fc 
Psaltis (1999) and Revnivtsev, Gilfanov, & Churazov (2001), 
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although those authors parameterized hardness by the pho- 
ton index, F, and they did not include HEXTE data in their 
fits. Observations in the 'hysteretic region', which have rel- 
atively similar frequencies, deviate from this trend slightly. 
The overall behaviour, however, is roughly consistent with 
the PSD frequencies being proportional to i^ z ^ (although 
the highest frequency exhibits a sl ightly flatter trend). We 
comment upon this further in §5.4 below. 



coronal compactness from the kotelp+gauss fits of Ta- 
ble AJ. Again, the plateau of the observations in the hys- 



We have also calculated the average time lag between 
the variability in the 0-5 keV and the 10-20 keV energy 
bands. In general, the broad band variability of GBHC 
in their hard state shows the hard variability lagging be- 
hind the soft variability, with the lag increasing both for 
lower Fourier frequencies and for larger energy separations 
(Miyamoto & Kitamoto 1989; Miyamoto et al. 1992; Nowak 
et al. 1999; Nowak, Wilms, fe Dove 1999; and references 
therein). As seen in Figs, [j^, [y| and [l2| the variability 
time lags span a large dynamic range from w 10 -3 sec at 
high Fourier frequencies, to « 0.1 sec at low frequency. Var- 
ious models have attributed this broad range of observed 
t ime delays to photo n diffusion in a very extended corona 



(Kazanas et al. 1997), wave propagation in a disc (Nowak 



et al. 1999c), or temporally correlated fl ares covering an ac- 
cretion disc ( Poutanen fc Fabian 1999| ) . In order to both 
increase the signal to noise and to describe the overall char- 
acteristic behaviour of the time delays, we have averaged the 
calculated time lags over a range of a factor of 5 in Fourier 
frequency. The resulting average time lags are plotted vs. 
3-9 keV flux in Fig. 0. 

At the highest frequency, we see that the time delays 
are roughly consistent with being at a common minimum 
value of ~ 10 -3 Hz. We have previously argued that such 
a minimum value of the time delay sets a lower l imit of 



roug hl y 30 G M/c for the size scale of any corona (|Nowak| 



et al. 1399c). For lower Fourier frequencies, we see that in 
the 'hysteretic region' time lags are relatively uniform from 
observation to observation. For observation P40108_03, co- 
incident with the 1999 soft to hard state transition, there 
is an increase in the characteristic time lag. A very similar, 
but more dramatic, time lag increase also has been seen in 
the soft /har d state transitions (or 'fa iled' state transitions) 
of Cyg X-l ([Pottschmidt et al. 2000|). At lower fluxes, the 



time lag is seen to drop. 

Nowak, Wilms, & Dove ( [l999b| ) have previously argued 
that the simultaneous decrease in the characteristic PSD fre- 
quencies and in the time lags seen for observation P20181_05 
(i.e., the lowest flux diamond points in the fi gures) argues 
against the flare models of Poutanen & Fabian (19991), where 
one expects the PSD time scales (i.e., inverse frequencies) 
and the time lags to be correlated. The trend exhibited by 
observation P20181_05 for the PSD time scales and time 
lags to be anti- correlated is seen to continue towards lower 
flux. Interestingly, the lowest flux observation shows that 
the power spectra averaged about 0.1 Hz and 4 Hz exhibit 
a lag of the soft variability behind the hard varia bility. We 
note, however, that as discussed by Nowak (200C), the net 



overall time lags could in fact be made up of several indepen- 
dent processes — which we would be averaging over in these 
figures — some of which have soft variability lagging hard 
variability, and some of which have hard variability lagging 
soft variability. 

In Fig. [ll] we show the time lags plotted against the 



teretic region are apparent. The overall trend is for less com- 
pact/softer spectra to show longer time lags than more com- 
pact/harder spectra, albeit the faintest two observations de- 
viate from this trend. This deviation for the faintest two ob- 
servations is similar to that seen for the anti-correlation be- 
tween 'reflection fraction' and compactness (from the eqpair 
fits) seen in Fig. |S| Plotting time lag vs. reflection fraction 
from the eqpair fits, we see that there is indeed a correlation 
between the two, as shown in Fig. [l| where higher reflec- 
tion is seen to imply longer time lags of the hard variability 
behind the soft variability. If the reflection fraction is indeed 
truly measuring the relative amount of cold material being 
illuminated by the corona, then perhaps these data are argu- 
ing for low compactness/soft spectra being associated with 
a larger coronal size, which would yield longer time lags due 
to increased propagation or diffusion path lengths while si- 
multaneously illuminating a larger reg ion of the cold disc. 
We comment upon this further in §tj.4l. 



4 SOFT STATE OBSERVATIONS 

The soft state observations discussed here, P40108_01 and 
P40108_02, occurred shortly prior to a return to the hard 
state. A discussion of brighter soft state observations of 
GX 339—4, taken in early 1998 after the initial transition 
from the hard to soft state, can be found in Belloni et 



al. ( |1999D . Observations P40108.01 and P40108.02 are very 
characteristic of 'classic' BHC soft states (Nowak 1995, and 
references therein). Both observations exhibit rms variabil- 
ity 3% between fa 10~ -1 Hz, and are totally dominated 
by Poisson noise residuals at higher frequencies. 

In terms of their spectra, both observations show a 
strong soft excess at energies 8keV, with a V ~ 2 power 
law at higher energies; see Fig. [l3| Such observations are typ- 
ically well-described by a multi-temperature disc blackbody 
with maximum temperature of « 1 keV, plus an additional 
power law, typically with a photon index <; 2. For an exam- 
ple of similar soft state BHC spectra observed with RXTE, 
we refer to our recent work on LMC X-l and LMC X-3 
flNowak et al. 200l[ [Wilms et al. 200l[). For these three soft 



state BHC, a broad and strong Fe line seems to be required 
(~ IkeV width, « 200-9 00 eV equivalent width). As dis- 
cussed by Nowak et al. ( |2001 ), there is some worry that 
such a line could be an artifact due to the phenomenological 
fit components, namely a multi-temperature disc blackbody 
and a power law, crossing over one another at an energy 
~ 6keV. A broad line could be required to remove an arti- 
ficially created "inflection point". 

As with our LMC X-l and LMC X-3 spectra, to assess 
this possibility for the soft spectra of GX 339—4, we turn to 
the compps model, and choose the "slab geometry" (i.e., a 
planar disc and corona, with seed photon injection occurring 
at the disc-corona interface at the midplane). As this model 
self-consistently calculates the high energy tail from the seed 
photon distribution and the coronal parameters, there is less 
worry about creating an artificial Fe line. In addition to 
the coronal parameters, we also allow for reflection from an 
ionized slab, with a relativistically smeared reflection profile. 
(We choose the /3 — 10 option in the compps model, and fix 
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Figure 11. Same as in Fig. |ro[ except here plotted vs. compactness from the kotelp+gauss fits of Table Al. Left: Average frequency 
0.7Hz. Middle: Average frequency 3.7Hz. Right: Average frequency 18Hz. 



the inner and outer disc edge to be 6 GM/c 2 and 10 3 GM/c 2 , 
respectively. The disc Fe abundance is set to 4, and the disc 
inclination angle is frozen to 45 °.) To m odel the line, we 



use the laor model in XSPEC (Laor 1991). Again, we set 
the inclination angle to 45°, and we freeze the inner and 
outer disc edge at their maximum values. We allow the line 
energy, normalization, and the disc emissivity index to be 
free parameters. 

The 



Results for these fits are presented in Table A3 



fits are very similar to our prior results for LMC X-l and 
LMC X-3. Specifically, the disc inner edge temperatures are 
fcldisc ~ 540 eV, the fitted coronae are very hot with kT e « 
130 keV, and they are also optically thin, with T es « 0.3. In 
addition, the fitted Fe lines have energies w 6.8 keV (con- 
sistent with ionized iron), are very strong (600 to 900 eV 

equivalent widths^), and are very broad. Emissivity indices 

t Note that we do not quote error bars for the line equivalent 
widths in Table [A3] , as the error bars are complicated functions 
of the other fit parameters, not just of the line normaliaa 



zatio 



General ly, the equivalent widths deviate from the best fit val- 
ues by approximately — 100 eV to +500 eV, depending upon what 
other parameters are allowed to vary. 



are > 5, indicating lines that are extremely skewed towards 
the inner edge of the disc. The fitted reflection fractions, 
in contrast to the line equivalent widths, are small, with 
Q./2-k 0.2. In this respect, these observations are similar 
to P40108_06 and P40108_07. The cold reflector, consistent 
with the line energies, is ionized. 

The unfolded spectrum for observation P40108_01 is 
presented in Fig. [H| One question that arises, is the broad, 
large equivalent width line an artifact of our assumed model? 
For example, could the higher energy line with a sharp blue 
edge (due to the emissivity index being > 5), in reality be 
mimicking an unsmeared, neutral reflection edge at 7.1 keV? 
(The smeared, ionized edge is broader and extends to higher 
energies.) This would be expected if the reflection were due 
to a combination of an ionized, relativistic disc and cold, 
neutral matter much farther away from the compact object 
(i.e., a flared outer disc edge, the companion star, etc.). We 
have tried adding a neutral edge at 7.1 keV, and this does re- 
duce the Fe line equivalent width, but only by ~ 300 eV. As 



for our LMC X-l and L MC X-3 observations (Nowak et al 



2001; Wilms et al. 2001), a fairly strong, broad line is a ubiq- 



uitous feature of all the models that we have attempted to 
fit to the GX 339—4 soft state data. 
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Figure 13. Unfolded energy spectrum of GX 339-4 (P40108.01) with the compps Comptonization model fitted to the PCA+HEXTE 
data. 



Th 3 other aspect that is to be noted about Fig. 


15 is the 


et al. 1997b) and r es <- 


possibility of a hardening at energies above ~80keV. Due to 


flGierlMski et al. 1997 



the very low hard flux, this hardening is entirely consistent 
with the uncertainties of the background measurement (see 
§5.6). We note also that we have used the compps model to 



fit non-thermal Comptonization models, and obtain equally 
good, purely non-thermal fits. Thus, these data are unable 
to distinguish between a thermal and non-thermal electron 
distribution for the Comptonizing electrons. 



5 DISCUSSION 

5.1 Comparison to Previous Comptonization 
Models 

As discussed in §[Tl], the compps, eqpair, and kotelp mod- 
els endeavor to be physically motivated, realistic descrip- 
tions of Comptonization in compact object systems. The 
results that we have found here with these models are in 
many ways comparable to previous analyses of GX 339—4, 
as well as the spectrally similar source, Cyg X-l. For exam- 
ple, previous work with versions of the compps model have 
found r cs ~ 1 and kT ~ 100 keV for Cyg X-l (Poutanen 



1 and kT ~ 50keV for GX 339-4 
Our own previous analysis of joint 
PCA+HEXTE data of Cyg X-l u sing the kotelp m odel 
found r cs ~ 1 and kT ~ 90keV ( [Dove et al. 1998j ). On 
the other hand, we previously have found somewhat lower 
temperature (kT ~ 20-40 keV) and higher optical depth 
(t cs ~ 3) coronae for GX 339-4 flWilms et al. 1999[ ). How- 
ever, unlike the aforementioned studies which were broad- 
band (« 2- > 200 keV), the work of Wilms et al. (]l99S|) 



fit the low energy (PCA) and high energy (HEXTE) data 
separately. Recent work fitting broad-band data of Cyg X-l 
(Maccarone & Coppi, in prep.) with the eqpair model sug- 
gests similar coronal parameters as discussed above, and in- 
terestingly enough suggests at times a large equivalent width 
line (200 eV) with a low reflection fraction (Q/2n ~ 0.2), 
similar to what we have found for several of our observa- 
tions of GX 339—4. Thus, there is at least general consen- 
sus that "realistic" Compton coronae models with r cs ~ 1, 
kT ~ 100 keV, and low reflection fractions, Q/2n ^ 0.5, pro- 
vide reasonable descriptions of hard state black holes. (See 
§^ for references to recent work concerning soft state black 
holes, as well as Frontera et al. 2001, Gierliriski et al. 1999.) 

The above studies, however, for the most part consid- 
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ered either single observations, or a fairly small group of 
observations. Studies that have considered lar ge groups of 



1999 



obse rvations analyzed in a s y stematic fashion ( Wilms et al 
1999 ;_ jdziarski et al. 



Revnivtsev et al. 2001 ) have 
tended to focus on more phenemonological models (e.g., 
power- laws with reflection) . In addition, some of these stud- 
ies have ignored or lacked high energy data above ^ 30 keV 
(e.g., Revnivtsev, Gilfanov, & Churazov 2001, Zdziarski, 
Lubihski, & Smith 1999). Within these studies, there has 
been disagreement as to the interpretation of the results, 
even given the same data (see, for example, Wilms et al. 
1999 and Revnivtsev, Gilfanov, & Churazov 2001). For this 



reason, as well as the concerns outlined in S 2.2 , in the fol- 



lowing sections we examine our results more carefully. We 
begin by comparing our data to observations of the Crab 
nebula and pulsar. 

5.2 Crab Ratios 

In Fig. |l4| we show the normalized ratio for observations 
P20056_07 and P20181_05 with respect to the Crab neb- 
ula and pulsar observation of 1997 September 05, and for 
observations P40108_03 and P40108_06 with respect to the 
Crab nebula plus pulsar observation of 1999 February 24. 
Compared to direct fitting of spectral models, such ratio 
plots have the advantage that they are relatively insensi- 
tive to the uncertainties in th e calibration of the response 
matrix ( santangelo et al. 1995 ) , although they ar e less sen- 



sitive to detailed features than spectral fitting (Kreyken 



bohm et al. 1999 ). The Crab ratios were generated by first 
subtracting the background from both the Crab and the 
GX 339—4 observations. The Crab spectrum was then mul- 
tiplied by £ r c.ab-r 3 39 where rCrab = 2.18 (see Appendix) 
and where F339 = 1.73 (i.e., the photon index for observa- 
tion P20181_04). All spectra were then normalized to the 
same flux level at 10 keV. The Crab ratio was then obtained 
by dividing the GX 339—4 pulse height analyzer (pha) data 
by the modified Crab spectrum. 

Fig. |l4] shows most of the features that we have also 
found in the detailed spectral fitting and thus provides us 
with additional confidence that the general behaviour seen 
was not due to the response matrix uncertainty. Higher 
fluxes clearly yield softer spectra. Most of this variation 
is due to the soft X-rays, while the hard spectrum shape 
stays approximately constant. This is completely consistent 
with the results we previously obtained by modelling some 
of these data with a broken power-law model, where we 
found that the variation in the power-law index at softer 
X-rays was nearly twice as great than that at higher en- 
ergies (specifically, see Tabl e 5 of Wilms et al. 1999). We 
comment on this further in £ 5.3 below. 



The Crab ratios also highlight the broad feature present 
at ~ 6.4 keV, which in many models is consistent with a 
broad Fe line. In fact, fits with realistic Comptonization 
models always seem to require such a broad Fe line, as we 
first noted with fits to RXTE data of Cyg X-l (Dove et al. 
1998; who used a very early version of the PCA response ma- 
trix). With every subsequent revision of the PCA response 
matrix, realistic Comptonization models of both Cyg X-l 
and GX 339—4 data have continued to requ ire the pres- 
ence of such a broad Fe line (Wilms et al. 1999). Taking the 
orthodox interpretation of the feature present in the Crab 



ratios as a broad Fe line, these ratios provide us with mor e 
confidence in some of the unusual results obtained in §2.5 

First, observation P20181_05 does indeed have the 
smallest equivalent width line, as was determined by the 
fits. In addition, the line is very narrow, and is essentially 
consistent with the PCA resolution. Also shown in Fig. [l4| is 
the fact that the P40108 observations (i.e., the observations 
taken after the return to the hard state in 1999) indeed do 
have larger equivalent width, and for the most part broader, 
Fe lines. This is true even into quiescence, as shown by the 
Crab ratio for observation P40108_07. Note also that de- 
spite the fact that observation P40108_06 and P20181_05 
have similar flux levels, similar fitted compactnesses, and 
similar fitted reflection fractions, the Crab ratio shows that 
P40108_06 has the clearly stronger and broader Fe line (cf. 
Revnivtsev, Gilfanov, & Churazov 2001). Again, however, 
as discussed in §2.5, the fitted physical line widths, <r, are 



complicated by one's assumptions for the line energy, the 
assumed reflection properties, and the continuum spectrum. 

The Crab division is too insensitive to reveal some of 
the finer effects found by our modelling. For example, one 
trend seen in the eqpair fits to the data is for the temper- 
ature of the seed photons to decrease with decreasing flux. 
In the simple sphere+disc coronal model, one expects the 
approximate relationship fcTbb o c ^ < 7 1//2 -F^, 4 , where Fbol is 



the total flux of the system (see §5.4 below). Thus, utilizing 
the results of Table [l] and Al , one expects an approximately 
factor of 2 change in the seed photon temperature from the 
brightest to the faintest observation, which is slightly larger 
than actually found for the P20181 and P20056 observations. 
Dramatic changes, however, are found for the first few ob- 
servations after the return to the hard state in 1999, as we 
show in Fig. [H]. The seed photon temperature is seen to be 
greatly lowered, and then slowly recovers to the previously 
observed trend over the course of the next several months of 
decline into quiescence. As discussed in such a low seed 
photon temperature was required to reproduce a soft excess 
at energies ^ 4keV, where we chose very small systematic 
errors, thus making the trend extremely statistically signifi- 
cant. As seen from Fig. |l4| however, fractionally this trend 
is subtle and small, and therefore we cannot use the Crab 
ratios to rule out the possibility of a systematic trend in the 
(time-dependent) PCA response matrices that we use. 

One argument in favor of these features not being sys- 
tematic is that the transition to the hard state coincides 
with the turn on of the radio emission, which initially ap- 
pears in observation P40108_03 with a mildly optically thin 
radio spectrum and large amplitude (tens of percent) v ari- 
ations on J; 10 minute time scales flCorbel et al. 2000j ). It 
is tempting to associate the unusual soft excess/low fitted 
seed photon temperatures with the initial formation of the 
radio outflow. 



5.3 Is the T-n/27r Correlation' Real? 

As discussed above, the Crab ratios clearly show that the 
GX 339—4 hard state spectra can be well-described by a bro- 
ken power law plus a broad line, with the low energy (softer) 
pow er law varying mor e than the high energy (harder) power 
law ( |Wilms et al. 1999| ). Using more physical models, this ef- 
fect can be found again in the relation between the reflection 



fraction and the slope of the underlying power law ( Zdziarski 
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Figure 14. The normalized Crab ratio for observations P20056_07 and P20181_05 with respect to the Crab observation of 1997 September 
5, and for observations P40108-03 and P40108-06 with respect to the Crab observation of 1999 February 24. See text for a complete 
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et al. 1998). Since the shape of the spectrum above lOkeV 
changes only slightly while the soft spectral index shows a 
strong variation, the comparably constant spectrum above 
10 keV is reproduced by having a stronger reflection compo- 
nent for brighter/softer spectra. This ability of a reflected 
component to phenomenologically ha rden a soft power law 
has caused some, including ourselves flWilms et al. 1999| ), to 
question the reality of the putative F-fl/2-K correlation. The 
Crab ratios, however, remove any doubt. Viewed from the 
purely phenomenological perspective of a broken power law 
model, the V-Q./2n correlation' is real^ 

Viewed from a physical perspective, the exact nature 
of these two power law components is still uncertain. As 
measured with the Oriented Scintillation Spectrometer Ex- 
periment (OSSE) on board the Compton Gamma Ray Ob- 
servatory (CGRO), the high energy spectrum of GX 339—4 
is very well fit by an exponentially cut-o ff power law, with - 
out any additional reflected component ( prove et al. 199*8 ). 
HEX TE spectra fitted b y themselves yield very similar re- 
sults (Wilms et al. 1999). Reflection models are not detect- 



ing any characteristic "spectral curvature" in the 30-100 keV 
band. Instead, the reflection models are primarily sensing 
the break between the low and high energy power laws, with 
a larger break requiring a greater reflection fraction. 

Even here, the exact value of the "reflection fraction" is 
su bject to much unc ertainty. As we first discussed for Cyg X- 
1 ( Dove et al. 1998 ), since its i=S lOkeV spectra can be well- 
described by a spectrum without reflection, the fitted reflec- 
tion fraction is strongly dependent upon one's assumptions 



3 Note that the soft state PC A spectra presented in §|4| can also be 
reasonably well r eprese nted by the same models used in Table 5 
of Wilms et al. (1999); i.e., a weak disc blackbody with peak 
temperature 



about the nature of the ^ 10 keV spectrum. For Cyg X-l, we 
saw that if we allowed a large equivalent width and broad 
Fe line, and a higher effective temperature than 'usual' soft 
excess, no reflection was required (again, using earlier ver- 
sions of the PC A responses matrices). On the other hand, 
the kotelp model, which does have an effective reflection 
fraction of Q./2ir ~ 0.3, also fit those data. This is similar 
to t he situation here, where the fits of GX 339—4 discussed 
in £2.3 — which allowed for high temperature seed photons 
and large equivalent width, broad lines — yielded syst em- 
atically different "reflection fractions" than the fits of §2.5 



250 eV, a gaussian line, and a broken power law. 
For those observations, however, the Fe line has an equivalent 
width « 1.5 keV, and the broken power law indices are ft: 4.4, 2.2 
for the lower and higher energy power laws, respectively. 



The same, physically realistic Comptonization model, which 
self-consistently calculates the Comptonized spectrum from 
the seed photon spectrum and coronal parameters, yields 
different, nearly equally as good, answers for different sets 
of starting values. Thus, it is very difficult to uniquely de- 
termine the 'best' continuum model and thereby obtain the 
most 'accurate' value for the reflection fraction (cf. di Salvo 
et al. 2001). 

Again, however, viewed purely from the phenomenolog- 
ical perspective of a broken power law, the relative trends 
are giving us information concerning the degree of the break 
between the low and high energy power laws. As shown in 
Fig. jl4| those observations that show little or no break in 
the power law also yield fitted reflection fractions close to 
zero. This is even given the fact that the Fe line, which may 
or may not be related to the same component that causes 
the 'reflection', does not vanish. 

Fig. [l4| also provides an explanation for one of the other 
results discu ssed by Wilms et al. (1999). Specifically, Wilms 
et al. ( |l999| ) showed that the P20181 observations discussed 
here could be fit by an ionized reflector model with a weak 
and narrow Fe line. The r-f2/27r correlation was absent, but 
a hardness-ionization parameter correlation was instead de- 
tected: softer spectra showed a more highly ionized reflector. 
Fitting a continuum model consisting of a soft excess (mod- 
elled as a 250 eV maximum temperature disc blackbody) , 
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Figure 15. Top: Seed photon temperature vs. 3-9 keV PCA flux, for the eqpair models. Bottom: Seed photon temperature vs. coronal 
compactness, for the eqpair models. 



and a power law fitting predominantly the 4~7keV data, the 
power law break can then be viewed as an edge, with the 
position of the edge shifting to higher energies (as would be 
appropriate for a more ionized reflector) for softer spectra, 
as is apparent in Fig. Such models quite acceptably de- 
scribe the data, and in fact can also be shown to describe the 
P20056 observations as well. Again, we see that given the 
poor spectral resolution of RXTE coupled with the inherent 
degeneracies of the spectral models (themselves consisting 
of broad, overlapping components), it is difficult to obtain a 
unique interpretation of the observed features. 



5.4 A (Too?) Simple Coronal Model 

In §131, we found that for the kotelp model fits the character- 

*i 3/2 

istic PSD frequencies scaled approximately oc £ c ■ Here we 
provide a simple hypothesis for this behaviour. The kotelp 
model fits a relative compactness of the corona as compared 
to that of the disc. Let us assume that the transition ra- 



dius between the disc and the corona occurs at a radius of 
Rt, while the bulk of the energy release from the corona is 
due to falling through the gravitational potential drop to 
the marginally stable orbit at R m . The compactness of the 
corona is proportional to the energy release in the corona 
divided by the coronal radius, i.e. oc _R m i? t _1 . The intrin- 
sic compactness of the disc, on the other hand, goes as its 
gravitational energy release (oc R^ 1 ) divided by its charac- 
teristic radius (Rt) and thus scales as oc R^~ 2 . The relative 
compactness of corona to disc is therefore Z c oc Rt/Rm- This 
is essentially what we fit with the kotelp model. 



The fits of 



showed that l c ranged from ~ 1 (for 
the observation after GX 339—4 returned to the hard state 
in 1999), to ~ 7 for some of the faintest observations. In 
the above simple model, the flux-hardness anti-correlation 
then can be ascribed to the corona-disc transition radius 
varying from R t » R m to R t « 7 R m <S 40 GM/c 2 , while 
the accretion rate drops by a factor of i; 4 (assuming that 
the 3-200 keV flux is a fair measure of the bolometric flux) . 
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If we take the view point of Psaltis & Norman (200C) 
that the observed PSD is made up of resonant time scales 
from the corona-disc transition region, many of these char- 
acteristic frequencies scale as R^ 3 ^ 2 and therefore should 
scale as Z<7 3 ^ 2 in the above scenario, consistent with the re- 
sults shown in Fig. ^| One implicit assumption made here is 
that the accretion rate, M , scales out of the ratio between 
the coronal and disc compactness. That is, we are assum- 
ing an efficient accretion disc. A further assumption is that 



the scaling of the frequencies is indeed oc i?, 



-3/2 



Whereas 



this is true for the dynamical, thermal (for fixed a), vertical 
epicyclic, and radial epicyclic (for R ^> 8 GM/c 2 ) frequen- 
cies, this is not true for the Lense-Thirring precession or 
viscous damping frequencies. 

Although the distance, absolute luminosity, and mass of 
GX 339—4 are not well-known (see Zdziarski et al. 1998 for 
a thorough discussion of these points), it has L 5% LEdd 
if its fractional Eddington luminosity is comparable to other 
hard state GBHC (Nowak 1995, and references therein). 
Thus, in any thin disc Shakura-Sunyaev model, whether it is 
gas or radiation pressure-dominated, one expects the viscous 
time scale to be several thousand seconds at R ~ 20 GM/c 2 , 
even for a relatively large. This is substantially longer than 
the low-frequency cutoff time scales presented in Fig. ^. In 
addition, the ratio between the viscous time scale and the 
thermal time scale is approximately oc (L/LEdd)~ 2 400, 
which is a larger dynamic range than that presented in 
Fig. ^. Let us assume, however, that the transition region 
between a corona and disc is fixed and is relatively a large 
fraction of the transition radius, on the order of H/R ~ 0.1- 
0.3. The dynamic range of frequencies shown in Fig. ^ then 
can be consistent with including the viscous time scale, and 
furthermore the viscous frequency would also scale as R^/ 3 ^ 2 . 

Within this simple model, compactness changes are 
achieved solely by varying the coronal radius, and not by 
changing the underlying geometry. This is in contr ast to 
the suggestion of Zdziarski, Lubihski, & Smith ( 1999] ), who 
hypothesize that the hardness and reflection fraction are 
regulated by the degree of 'overlap' between the cold disc 
and quasi-spherical corona (more overlap yields both softer 
spectra an d greater reflectio n fraction). We had previously 
suggested (Dove et al. 1998) that some amount of overlap 



was required to produce the broad, large equivalent width 
line that was required for kotelp models of both Cyg X-l 
and GX 339—4. We have attempted such simulations with 
the kotelp code; however, our results indicate that the coro- 
nal region overlapping the disc cools to the inverse Compton 
temperature, kTic ~ 4keV, and thus fails to yield the re- 
quired additional reflection. This does not mean such mod- 
els are ruled out, rather that one needs to carefully model 
additional, substantial heating in the overlap region in or- 
der to maintain the local coronal temperature. It has been 
suggested that such local heating can ar ise from thermal 



conduc tion betw een the corona and disc (Meyer & Meyer- 



Hofmei iter 1994), although it has not been rigorously shown 
that such heating is sufficient. 

The above 'overlap' model has been suggested as possi- 
bly being applicable to an overall T-fi/2-Tr correlation' com- 
mon to Seyfert 1 galaxi es and galactic black hole candidates 



(Zdziarski et al. 1999). However, as discussed above, the 
viscous time scale in a thin, low luminosity accretion disc 
about a low mass black hole is hundreds to thousands of 



seconds. Thus, each individual observation discussed here 
is integrated over many viscous time scales, while different 
observations can, and likely do, represent vastly different 
intrinsic accretion rates through the disc. This is in con- 
trast to a 10 s M Q black hole in a Seyfert 1 (e.g., NGC 5548; 
see Chiang et al. 2000, and references therein), which likely 
has a characteristic viscous time scale of years. The week to 
week spectr al variations of NGC 5548 discussed by Chiang 
et al. (2001) therefore correspond to thermal time scales and 
faster. Such thermal time scales are integrated over in a few 
seconds in our observations of GX 339—4. Thus, whereas 
their are undoubtedly commonalities between Seyfert l's 
and GBHC, it would be surprising if the exact same models 
were responsible for the spectral correlations seen in both. 

A major failure of the above simple model is that it 
does not explain the putative reflection fraction correlation 
with compactness. To consider this further, we turn to the 
eqpair model results, which do allow for reflection changes. 
The eqpair model fits agree in rough outline, but not in 
detail, with the kotelp model fits. Specifically, both require 
coronal temperatures 150 keV (although eqpair fits tend 
to have optical depths lower by a factor of approximately 2- 
3), and both are consistent with overall low reflection frac- 
tions of £l/2n ^ 0.5. This reflection fraction, however, is 
clearly variable in the eqpair fits. In addition, the eqpair 
model also requires a seed photon temperature variation, 
which in turn relates to the "effective coronal radius" for 
these models. Since we have fixed the seed photon compact- 
ness to one, the fitted eqpair coronal compactness scales as 
Fboi/(o'SBT^ h R 2 ), where F bo i is the bolometric flux, ass is 
the Stefan-Boltzman constant, Tbb is the temperature of the 
seed photons, and R-, is the "effective radius" for the seed 
photon distribution, which for the eqpair model is assumed 
to be on the order of the coronal radius. Assuming that the 
3-200 keV flux is a fair measure of the bolometric flux, and 
taking a distance of 4 kpc and a mass of 4 Mq for GX 339—4, 
this effective radius, in units of GM/c 2 , is approximately 
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Note that in the hysteretic/overlap region (the bright- 
est seven P20181 observations), the effective radius and fre- 
quencies are roughly constant and, within the large uncer- 
tainties, ^ 100 GM/c 2 . However, for the 1999 return to the 
hard state, the radius trends are exactly opposite of those 
implied by the kotelp fits. Specifically, the effective radius 
is seen to be very large, and then shrinks. This behaviour, 
although counter to an association of the PSD frequencies 
with a coronal radius, is in agreement with the association 
of the time laps with a characteristic size scale, as we also 
show in Fig. [16| Is this necessarily in contradiction to the 
kotelp results? 

For the 1999 hard state observations, the large reduced 
X 2 's achieved for most of the kotelp models were primarily 
due to residuals in the 3-4 keV band (with some contribu- 
tion from an "edge- like" residual near lOkeV). Otherwise, 
the kotelp models fit the data reasonably well, especially at 
high energies. The low seed photon temperatures found with 
the eqpair model were largely driven by the 3-4 keV resid- 
uals. We speculate here that the "true" model might be a 
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Figure 16. Top: Characteristic PSD frequencies vs. "effective radius" for the seed photon distribution from the eqpair fits. Bottom: 
Time lags between soft and hard variability (average PSD frequency 0.7 Hz) vs. effective radius from the eqpair fits. For clarity, only 
the higher signal-to-noise P20181 and P40108 observations are shown. 



combination of the suggested compactness/frequency /time- 
lag/radii relations above. That is, perhaps the correct model 
involves the radii of the base of the corona increasing with 
decreasing flux/increasing compactness, while the vertical 
extent of the corona decreases for decreasing flux/increasing 
compactness. The initial transition to the hard state, exhib- 
ited by observation P40108_03, could be associated with the 
formation of a large vertical scale jet, with X-ray emission 
extending to nearly 1000 GM/c 2 , while its base has a radius 
on the order of the marginally stable orbit. The character- 
istic PSD frequencies could then be generated within the 
disc, while the characteristic time lags could be generated 
via propagation along the length of the jet. An extended 
jet would also "view" a much larger extent of the disc, and 
thereby have a greater reflection fraction. (The slightly low 
reflection fraction exhibited by observation P40108_03 given 
its long time lag, as seen in Fig. [E| might also be attributed 
to some amount of beaming away from the disc, as sug- 



gested, e.g., by Beloborodov 1999.) The reflection fraction 
would naturally decrease as the jet height decreases and the 
jet base widens. 

Clearly the above suggestion requires detailed theoret- 
ical modelling to verify. Furthermore, it is subject to great 
observational uncertainty. Specific ally , are the seed photon 
temperature trends shown in Fig. real or systematic er- 
rors? However, we note that the large effective radius implied 
by the eqpair fit to observation P40108_03 implies the pos- 
sibility for this observation, and the subsequent two observa- 
tions, of having a large contribution from synchrotron seed 
photons, as opposed to solely thermal seed photons. Taking 
the fitted eqpair coronal temperature and optical depths, 
and hypothesizing a magnetic field in equipartition with the 
thermal energy of the corona, then the expected flux (at 
4 kpc, assuming a 4 Mq black hole) due to synchrotron emis- 
sion within the corona is on the order of 
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an implied "seed photon" flux of ~ 1CP 9 
see that the synchrotron flux only becomes comparable to 
the thermal seed photon flux for an effective coronal radius 
ii 300 GM/c 2 , i.e. the implied effective radius of observa- 
tion P40108_03. Again, this observation, coincident with the 
initial transition to the hard state, showed unusually large 
radio variability on 10 minute time scales, a s well as ex- 



hibi ted a quasi-optically thin radio spectrum (Corbeletal 
200c()~Jhus, the low fitted seed photon temperature for the 
eqpair model may have been indicating the need to consider 
synchrotron seed photons. 



5.5 Implications for ADAFs 

An Advection Dominated Accretion Flow was first suggested 
as a m odel for the hard state spectrum of Cyg X-l by Ichi- 
maru (1977). Since that time, there has been much further 



research and applications to spectra of GBHC. One of the 
key questions in applying ADAF models has been, where 
is the transition between the hot, geometrically thick advec- 
tive flow and the cool, geomet rically thin outer disc located? 
Esin, McClintock, & Narayan ( 1997 ) had originally proposed 
that the soft to hard state change in GBHC was associated 
with an increase of the "transi tion r adius" from ~ 10 GM/c 2 
to w 10 4 GM/c 2 . Esin et al. ( 1998 ) later revised this range, 
for models of Cyg X-l, to » 60 -* 400 GM/c 2 . Recent mod- 
els of XTE J1 118+480 require the transition radius to be 
> 100 GM/c 2 flEsin et al. 200l| ). 

This aspect of the ADAF model can be tested with these 
observations. We see that for both the kotelp and eqpair 
models, with the exception of the first few observations af- 
ter the 1999 return to the hard state, the implied coronal 
radii are all <> 100 GM/c 2 . There is no evidence for large 
amplitu de variations in this radiu s as the observed flux de- 
creases (di Matteo & Psaltis 1999). In fact, both the eqpair 
and kotelp models suggest that for the faintest observation, 
P40108_07, which is 15 times fainter than the brightest hard 
state and 5 times fainter than the initial transition back to 
the hard state in 1999, that the effective coronal radius has 
actually shrunk in comparison to the previous hard state 
observations at higher fluxes. 

Evidence for a small coronal radius for observation 
P40108_07 also comes from the observed Fe line, which in 
Fig. |l4| clearly exhibits a broad, red tail. Fitting the 3-20 keV 
PC A data of this observation with a co mbination of a powe r 
law and a relativistic diskline model ( Fabian et al. 1989 ), 
we find the following limits on the inner radius of this emis- 
sion. Assuming a line at 6.4 keV, with a 45° inclination for 
the disc, and a line emissivity that is oc R~ 3 , the best fit 
inner radius is 14±| 3 GM/c 2 (90% confidence level). Re- 
placing the power law with a power law plus reflection (i.e., 
Magdziarz & Zdziarski 1995; here with a best fit value of 
n/2-K = 0.2), we find an inner disc radius of 6± 6 GM/c 2 . 
Replacing the reflected power law with an eqpair model as 
in 5E.EI, the inner disc radius is 13+f GM/c 2 . Such small 
implied inner disc radii are strongly counter to what would 
be the best evidence for the "typical" ADAF model. 

Of course, it is possible that the advection dominated 



.region only exists at radii Ss 100 GM/c 2 . The arguments pre- 
sented above about the scaling of the PSD time scales with 
the kotelp compactness parameter, which assumed efficient 
accretion, are only very weak, indirect evidence against ad- 
vection domination. However, none of the observations here 
require advection domination either, and there is no com- 
pelling evidence for any efficiency change from the soft to 
hard state spectra. The observations are all consistent with 
the spectral changes solely being determined by accretion 
rate changes. (Further evidence of this comes from the fact 
that the radio flux closely linearly tracks the X-ray flux, all 
the way into quiescence; Corbel et al. 2001, in preparation.) 



5.6 Alternative Models? 

As shown in Fig. some of our observations (most no- 
tably P20181_01) exhibit a hardening above ~100keV in 
the HEXTE data. One possible explanation would be that 
the hardening is hinting at presence of a high-energy power - 
law tail as seen, e.g., in Cyg X-l (McConnell et al. 2000). 



Power-law tails are generally attributed to the presence 
of a non -thermal electron population in the Comptonizing 



plasma (Bednarek et al. 1990 



Gierliriski et al. 1997; Copp 



199S). Alternatively, weak power- law tails have also been ex- 
plained in terms of thermal Compton ization in a Maxwellia n 
plasma with a temperature gradient ([Skib" o fe Dermer 1995|) . 
In these models, the temperature gradient serves to make 
the "effective" electron distribution non-thermal. See Mc- 



Connell et al. (200C) for a discussion of these mechanisms. 
In galactic black holes, power-law tails are generally associ- 
ated with the soft state or with transitions into the soft state. 
As these transiti ons are typic ally associated with enhanced 
radio emission ( |Fcnder 20o"o| ), one is tempted to see the 
presence of the power-law tail as an indication of the non- 
thermal el ectron distribution in a radio outflow (Markofi 
et al. 200l| ). 

A possible alternative explanation for the power-law 
tail, however, is that it is due to instrumental effects: The 
tail could also be due to an incomplete background subtrac- 
tion that might be caused either by the uncertainty in the 
estimation of the HEXTE live-time or by a transient back- 
ground source. Such checks are especially important since 
the HEXTE background contains two strong line complexes 
below 100 keV. In order to test for the presence of a back- 
ground source we extracted individual HEXTE background 
spectra for each of the two background dwells of the HEXTE 
clusters and then compared these spectra to each other. In 
the case of observation P20181_01, the difference between 
the plus and minus background position for cluster B is 
3.6 x 10~ 2 cps, consistent with statistical fluctuations. On 
the other hand, for cluster A the difference between the two 
background dwells amounts to almost 1 cps. The difference 
spectrum between the two dwells, however, does not resem- 
ble the observed hardening in GX 339—4. Furthermore, we 
note that we see the tail in both HEXTE clusters which 
makes it quite improbable that this 1 cps difference is the 
cause of the hard tail. We conclude that the question of the 
hard tail in GX 339—4 is still open for discussion, but that 
its association with a problem of the background subtraction 
is not strongly suggested. Observations with more sensitive 
gamma-ray detectors are needed to determine whether this 
tail is a real feature. 
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6 SUMMARY 

We have presented spectral and timing analysis of 22 in- 
dividual observations of the galactic black hole candidate 
GX 339—4, as observed by the Rossi X-ray Timing Explorer. 
Several of these observations were coincident with Advanced 
Satellite for Cosmology and Astrophysics, radio, and/or op- 
tical observations. Our chief results for the X-ray analyses 
are as follows. 

• As has been previously observed, the soft state has a 
strong soft excess below ~ 7keV, a weak hard tail, and 
very little X-ray variability. Also in agreement with previous 
observations, the hard state, although exhibiting a spectral 
break at w lOkeV, does not show a strong soft excess, does 
have a strong hard energy tail, and exhibits large amplitude, 
rapid X-ray variability. 

• Generally, the hard state exists at lower bolometric lu- 
minosities than the soft state; however, in terms of 3-9 keV 
flux, there appears to be a " hysteresis" in the state tran- 
sitions, as first discussed by (Miyamoto et al. 1995). The 



degree of this hysteresis in terms of the bolometric flux is 
unclear. 

• Brighter flux states tend to show softer overall spectra; 
however, the first few observations after the 1999 return to 
the hard state are unusually soft. (This is the first "return" 
to a hard state observed by RXTE in GX 339-4.) In addi- 
tion, we have discovered that the trend of hardening spectra 
with decreasing flux is broken for the lowest flux observation. 

• Both the soft state and hard state spectra can be rea- 
sonably well represented by Comptonization models, with 
parameters similar to previous analyses. However, there is a 
great deal of degeneracy in these models, with different mod- 
els producing comparable fits. The kotelp and eqpair hard 
state models both require coronal temperatures <; 100 keV 
and both are consistent with reflection fractions Q./2-n 0.5; 
however, the former requires an approximately factor of 
three greater optical depth. In addition, there are basic dif- 
ferences between these models (e.g., 100% of the seed pho- 



tons pass through the corona in the eqpair model, whereas 
~ 30% of the seed photons pass through the corona in the 
kotelp model). 

• The hard state spectra show a break at » 10 keV, with 
the spectrum hardening at higher energies. Modelling this 
as a broken power law, the degree of this break increases 
with brighter/softer spectra. This has been dubbed the T- 
Q./2-k correlation' in previous analyses with reflection mod- 
els. Within the context of reflection models, in contrast to 
previous analyses, we find that the "reflection fraction" is 
better correlated with 3-9 keV flux than with the spectral 
hardness. 

• Most models require the presence of a broad, large 
equivalent width Fe line. For the 1997 hard state, the equiv- 
alent width of this line is correlated with the reflection frac- 
tion. For the 1999 hard state, unlike for previous obser- 
vations, the line remains broad and at relatively constant 
equivalent width, even into quiescence. This is counter to 
the simplest expectations of ADAF models. Ratios of the 
GX 339—4 data to observations of the Crab nebula and pul- 
sar show that these results represent real effects, and not 
systematic errors. 

• Timing analysis shows that brighter, softer spectra tend 
to have higher characteristic PSD frequencies than fainter, 
harder spectra, as has been previously observed in both 
GX 339—4 and Cyg X-l. We have, however, included more 
variability components than previous analyses, and further- 
more we have correlated our results with coronal parameters, 
rather than a power-law index. We here have discovered also 
that the variability time lags show the opposite trends from 
the PSD frequencies. Bright, soft spectra tend to show the 
longest time lag of the hard X-ray variability behind the soft 
X-ray variability. 

• A correlation between time lag and 'reflection fraction' 
was also newly discovered. The longest time lag was seen 
for the first observation taken after the 1999 return to the 
hard state, similar to observations of 'failed state transitions' 



in Cyg X-l (Pottschmidt et al. 200C). The 1999 return of 
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GX 339—4 to its hard state also exhibited unusual radio 
properties. 

• Variability properties, both for the characteristic PSD 
frequencies and measured time lags, were most uniform in 
the "hysteretic" range of 3-9 keV fluxes (cf. Revnivtsev, Gil- 
fanov, & Churazov 2001). 

• We have hypothesized that a radiatively efficient 
corona, wherein the base of the corona/jet decreases and 
its height increases, with increasing flux/softness, might be 
able to explain a number of the observed spectral-variability 
correlations. 

• More speculatively (due to uncertainties in instrument 
backgrounds and calibration), there is new evidence within 
the brightest, softest observations for a power law extending 
beyond the spectral roll-over usually associated with ther- 
mal Comptonization. The first few observations after the 
1999 return to the hard state also exhibited an unusual soft 
excess, which we speculated might be related to the turn on 
of the radio jet. 

Our own inclination is to believe that the hard tail observed 
for observation P20181_01 and that the soft excesses ob- 
served for observations P40108_03-_05 represent real phe- 
nomena, and not systematic effects. Their full implications 
can only be addressed with improved observations. For the 
former, there is some hope that Integral will provide re- 
fined high energy observations. For the latter, simultane- 
ous XMM-Newton/RXTE or Chandra/RXTE observations 
will be crucial for describing the seed photon regime of 
the spectra. This will be especially important, for exam- 
ple, in breaking the degeneracy between the 600 eV seed 
photon/50 keV corona fits (§[2.3 j) and the 100 eV seed pho- 
ton/150 keV corona fits (§ |2.5| ) given by the eqpair models. 
Furthermore, the ratio of the Comptonized flux to the seed 
photon flux (e.g., the "amplification factor") has implica- 
tions for the coronal geometry. In this sense, the eqpair 
and kotelp models are very different, and in fact, the op- 
tical (seed photon) to X-ray (Comptonized) flux might be 
constraining Seyfert 1 geometries in this manner (Chiang & 
Blaes 2001; in preparation). Finally, none of the currently 
available Comptonization models are truly adequate for de- 
scribing these data. Theoretical work, accounting for jet-like 
geometries and the presence of the known radio emitting 
outflow need to be more fully explored. 
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APPENDIX A: DATA ANALYSIS METHODS 
Al RXTE Data Analysis 

All RXTE results in this paper were obtained using the stan- 
dard RXTE data analysis software, f tools version 5.0 with 



the corresponding PCA response matrix, and HEXTE re- 
sponse matrix dated 97mar20c. Data selection criteria were 
that the source elevation was larger than 10° above the earth 
limb and data measured within 30 minutes of passages of 
the South Atlantic Anomaly or during times of high parti- 
cle background (as expressed by the "electron ratio" being 
greater than 0.1) were ignored. For one archival observa- 
tion, P20056_04, these criteria produced only 45 sec of us- 
able data, in contrast to the discussion of th is ob servation 
found in Revnivtsev, Gilfanov, & Churazov (2001). Analy- 



sis of archival Crab observations indicated that low residual, 
uniform results (in terms of measured spectral slope of the 
Crab from observation to observation) were best achieved by 
restricting the analysis to the first anode layer of the pro- 
portional counter units (PCUs). We therefore only consider 
this top layer throughout this work. For both spectral and 
timing analyses, we combined the data from all active PCUs 
(typically four or five PCUs were active at any given time). 

For spectral fitting, we limited the energy range of the 
PCA data from 3 to 22keV and the energy range of the 
HEXTE from 19 to 198 keV. PCA channels were grouped 
so that each had a minimum of 20 counts, while HEXTE 
channels were grouped by: a factor of 3 for channels 16-20, 
a factor of 5 for channels 22-51, a factor of 15 for chan- 
nels 52-126, and a factor of 24 for channels 127-198. To 
take into account the calibration uncertainty of the PCA 
we applied channel dependent systematic uncertain ties i n 
a manner similar to that described by Wilms et al. (1999). 
These uncertainties were determined from a power-law fit to 
observations of the Crab nebula and pulsar. As described on 
the RXTE Guest Observer Facility web pages, the high volt- 
age gain of the PCA detectors was changed in mid-March 
1999 (demarking the change from PCA 'Gain Epoch 3' to 
'Gain Epoch 4'). For all the P20181 and P20056 observa- 
tions, as well as observations P40108_01-_04 ('Epoch 3'), we 
applied 0.3% systematic errors in PCA channels 0-10, 1% 
systematic errors in channels 11-18, 0.5% systematic errors 
in channels 19-29, 1% systematic errors in channels 30-51, 
and 2% systematic errors in channels 52-128. For observa- 
tions P40108_05-_07 ('Epoch 4'), these same systematic er- 
rors were applied instead to PCA channels 0-8. 9-15, 16-24, 
25-43, and 44-128, respectively. Background subtraction of 
the PCA data was performed using the 'SkyVLE' model, as 
for our previous studies of GX 339—4 (Wilms et al. 1999). 

As discussed by Wilms et al. (1999), there is a system- 
atic difference between the PCA and the HEXTE detectors 
for the best fit photon index for observations of the Crab 
nebula and pulsar. PCA fits yield a Crab photon index of 
T = 2.18, while HEXTE fits yield T = 2.08. We have made 
no attempt to account for any slope differences; however, 
we have accounted for the cross-calibration normalization 
uncertainties of the PCA and the HEXTE instruments rel- 
ative to each other by introducing a multiplicative constant 
for each detector in all of our fits. Typically, if the PCA con- 
stant was fixed at 1, the multiplicative constants were ~ 0.9 
for the HEXTE A cluster, and a few percent lower for the 
HEXTE B cluster. 

In all of our spectral fits, we have modelled the neu- 
tral hydrogen c olumn using the cross sections of Morrison 
& McCammon (1983). Furthermore, unless otherwise noted, 
we have fixed the value for GX 339-4 to 6 x 10~ 21 cm" 2 . 
In the parameter tables below, italicized parameter values 
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indicate that the parameter was frozen at that value for that 
analysis. 

For the v ariabili ty analyses, we follow the methods of 
Nowak et al. ( 1999a ), and references therein. Throughout 
this work we have adopted the PSD normalization of Belloni 
& Hasinger (199C), where integrating over positive Fourier 
frequencies yields the root mean square (rms) variability of 
the lightcurve. We created lightcurves with a time resolution 
of 2~ 8 sec; therefore, the Nyquist frequency was 128 Hz for 
all observations. For the PSD fits discussed above, however, 
we fit the data from 3 x 10 -3 Hz to an u pper cutoff frequency, 

for each observa- 



which is listed in Table A4 



and Table A 5 



tion. Furthermore, we logarithmically binned each PSD over 
a range of Fourier frequencies, A///, that varied between 
0.06-0.15, depending upon the statistics of each observation. 

We performed variability analysis between two energy 
bands: 0-5 keV and 10-20 keV. These energy bands cor- 
responded to PCA pha channels 0-13 and 27-53 for the 
P20181 observations and for observations P40108_03-_04, 
and to PCA pha channels 0-10 and 22-44 for observa- 
tions P40108_05-_07. The P20056 observations used a data 
mode with very limited high energy resolution, therefore the 
high energy band for those observations corresponds to 10- 
100 keV, i.e., PCA pha channels 26-249. Since this is such a 
broad range, which includes many noise dominated channels, 
we did not calculate the time lags for the P20056 observa- 
tions. 



A2 ASCA Data Extraction 

We extracted data from the two solid state detectors (SIS0, 
SIS1) and the two gas detectors (GIS2, GIS3) on board 
ASCA by using the standard ftools as descri bed in the 
ASCA Data Reduction Guide ( |Day et al. 1998[ ). We chose 
circular extraction regions with radii of ~ 4arcmin for the 
SIS detectors, and « 6arcmin for the GIS detectors. We 
excluded approximately the central 1 arcmin of the SIS de- 
tectors to avoid the possibility of photon pile-up. We used 
the sisclean and gisclean tools (with default values) to 
remove hot and flickering pixels. We filtered the dat a wit h 
the strict cleaning criteria outlined by Brandt et al. (1996); 
however, we took the larger value of 7GeV/c for the rigid- 
ity. The background was measured from rectangular regions 
on the two edges of the chip farthest from the source (SIS 
data), or from annuli with inner radii > 6 arcmin (GIS data). 
These data were cleaned and filtered in the same manner as 
the source files. 

We combined the two SIS detectors into a single spec- 
trum, and we combined the two GIS detectors into a sin- 
gle spectrum, properly weighting the response matrices and 
backgrounds. Furthermore, we rebinned the spectral files as 
follows. For the SIS spectra, the channels were grouped: by a 
factor of 5 for channels 17-246, by a factor of 10 for channels 
247-266, by a factor of 15 for channels 267-311, by a factor 
of 20 for channels 312-351, and by a factor of 29 for channels 
352-380. For the GIS spectra, the channels were grouped: 
by a factor of 10 for channels 18-617, by a factor of 15 for 
channels 618-752, by a factor of 30 for channels 753-932, 
and by a factor of 39 for channels 933-971. Note that these 
observations were taken rather late in the lifetime of ASCA, 
after the SIS detectors in particular had sustained heavy ra- 
diation damage. Although the SIS and GIS detectors were 



in good agreement with one another from 1.8-10 keV, they 
were noticeably different from one another in the 1-1.8 keV 
range. (The GIS spectra were systematically softer than the 
SIS spectra in this energy range.) Taking this as possibly 
due to SIS radiation damage, we only considered SIS data 
in the 2-10 keV band, whereas we considered GIS data in 
the 1-10 keV band. 

We accounted for the cross-calibration uncertainties of 
the SIS and GIS instruments relative to each other and rel- 
ative to RXTE by introducing a multiplicative constant for 
each detector in all of our fits. Note that, relative to ob- 
servations of the Crab nebula and pulsar, ASCA fluxes are 
systematically lower than RXTE fluxes by ft! 80%, which we 
have confirmed with our simultaneous ASCA /RXTE obser- 
vations. Throughout this work we always quote PCA fluxes; 
however, this needs to be borne in mind when comparing the 
flux levels discussed here t o the ASCA flux levels presented 
in Table 1 of Wilms et al. (|l999|). 
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Table Al. Parameters for models fit to RXTE (PCA+HEXTE) hard state data of GX 339-4. Left: kotelp; 72 degrees of freedom 
(DoF) for P20181.01-P40108.04; 66 DoF for P40108.05-.07. Right: eqpair; 70 DoF for P20181.01-P40108.04; 64 DoF for P40108.05-.07. 
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Table A2. Parameters for near face on kotelp plus dust scattering halo models fit to hard state data of GX 339—4. 
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GIS, SIS, PCA, HEXTE 



Table A3. Parameters for the compps model fit to RXTE (PCA+HEXTE) soft state data of GX 339-4. 



Obs. 


kT c 
(keV) 


7c s 


kT dbh 
(keV) 


U/2tt 

(10- 2 ) 


(10 4 ) 


Eline 

(keV) 


^linc 


EW 

(keV) 


x 2 /DoF 


P40108.01 
P40108.02 


1231^ 
136l 8 


o-3ol2:2 2 
o.28i2° 2 


54+0- 01 
u - J ^-o.oi 

54+0- 01 
u -° -0.01 


1911 
20l! 


o-8l2: 4 3 


6.8i° :i 
6.812-.1 


5 4+°' 7 

5 9+0.5 
°' Z -0.4 


0.6 
0.9 


59/67 
68/67 
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Table A4. Fits to the 10-20 keV PSD of GX 339-4 (1997 hard state). 



Obs. / max (Hz) zfc-Lon QPOi QPQ 2 QPQ 3 zfc-Lor 2 X 2 /DoF 



P20181.01 


20 


A 




0.13+°°; 

— 0.01 


R 




o.o6+°-°; 

— 0.01 


R 




0.15+°° n 

— 0.00 






A 




6.7+°' 2 

— 0.4 


X 


10" 


■3 


133/97 






f 

J 




0.14+°°i 

—0.02 


f 




0.33+°'°° 

—0.00 


f 

J 




0.43+°°° 

— 0.03 






f 

J 




3.7+°i 

J ' —0.2 




















Q 






Q 




1.9+°i 

— 0.2 




















P20181.02 


20 


A 




0.15+°°; 

— 0.01 


R 




0.06+°°° 

— 0.01 


R 




—0.00 


R -- 


= o.o5+°°; 

— 0.01 


A 




7.0+°'? 

— 0.2 


X 


10~ 


3 


113/94 






f 




0.15+°°; 

—0.01 


f 




0.31+°'°? 

—0.01 


f 




0.42+°°? 
^ — 0.01 


f = 


= 0.81+°'°« 

—0.05 


f 




3.5+°-; 

— 0. 1 




















Q 




17±l 


Q 




cr+0.3 

— 0.3 


Q -- 


= 3.s+;- n 

— 1.0 
















P20181.03 


15 


A 




0.15+°°; 

— 0.01 


R. 




0.04+°-°; 

— 0.01 


R. 




0.17+°" 

— 0.00 






A 




7.1+°'? 

— 0.2 


X 


10" 


3 


74/56 






f 




0.11+°°; 

—0.01 


f 




0.30+°'°° 

-0.00 


f 




0.35+°°; 

-0.01 






f 




3.4+°? 

—0.1 




















O 




21+* 







i.6+°-; 

—0.1 




















P20056_01 


15 


A 




09+°'°i 
u u -0.02 


R. 




07+°' 01 , 
u -0.02 


R 




10+°°; 

u — 0.03 






A 




4 4+O.8 
-0.3 


x 


10~ 


3 


37/52 






f 

J 




0.13+°°2 

u — 0.02 


f 




0.37+°'°? 

—0.01 


f 

J 




0.49+°'°? 

—0.05 






f 

J 




3.4+°? 

J —0.5 




















Q 




8.1+2'2 

— 2.7 


Q 




1.9+°'5 

-0.5 




















P20056.02 


15 


A 




o.o6+°°; 

w ^ — o.oi 


R. 




0.15+°°; 

v — 0.01 












A 




3.5+°1 

—0.3 


X 


10~ 


■3 


67/73 






f 

J 




0.14+°'°8 
—0.03 


f 




0.40+°'°S 

—0.05 












f 

J 




4.0+°'5 




















Q 




i.5+°i 

— 0.3 


























P20056.03 


15 


A 




0.07+°°; 

— 0.02 


R 




0.16+°-°; 

— 0.01 












A 




r+0-6 
—0.3 


X 


10~ 


3 


54/41 






f 




0.10+°°? 

—0.02 


f 




0.40+°'°? 

—0.03 












f 




4.2+°'4 
—0.4 




















Q 




1.5+°i 

— 0.2 


























P20056.05 


10 


A 




0.04+°° n 

w u -0.00 


R. 




o.i3+°-°; 

w °— 0.01 












A 




2.9+°'o 
^—0.2 


x 


10" 


3 


54/38 






f 




0.16+°°I 

—0.03 


f 




0.54+°'°2 

—0.04 












f 




4.5+°° 
°— 0.8 




















Q 




2.1±g"| 


























P20056.07 


15 


A 




0.05+°°; 

— 0.01 


R. 




0.04+°-°; 

— 0.01 


R 




0.11+nnn 
— 0.00 






A 




2 5+°' 3 


X 


10" 


■3 


58/70 






f 




0.17+°'° 2 

—0.04 


f 




0.57+°°? 

—0.02 


f 




0.65+ n '° 3 

— 0.05 






f 




4 7 +o.6 

^•'-0.5 




















Q 






Q 




2.0+°1 
— 0.3 




















P20056_08 


10 


A 




0.04+°°; 

— 0.01 


R. 




0.12+°°; 

— 0.01 












A 




9 1+0.3 


X 


10- 


3 


26/29 






f 

J 




0.19+°°o 

—0.02 


f 




0.67+°'°! 

—0.04 












f 

J 




5 1+ 1 ' 6 




















Q 




2.0±° l 


























P20181.04 


15 


A 




0.15+°°; 


R 


_ 


o.o5+°°; 


R 




0.14+g-gJ 






A 




6 l+°' 2 


X 


10- 


■3 


78/56 






f 


— 


0.14l°;°; 


f 




0.34l°;°; 


f 




0.42l°;° 2 






f 


— 


fi +0.2 
J.D_ 2 




















Q 




18+ 5 


Q 




1 7+0.2 

— 0.1 




















P20181.05 


8 


A 




0.55+°;; 


R 


— 


0.26+°°i 












A 




13.5+J'2 x 10 


-3 


50/36 






/ 


— 


o.o2i° ; ( ( ;; 


f 




0-06±°'g 












f 


— 


9 1+0.2 




















Q 




6+ 01 
u - -0.0 


























P20181.06 


30 


A 




0.181°;°; 


R 




n (1^+0.02 


R 




17+ - 01 
u - ± '-0.01 






A 




7 0+0.5 


X 


10- 


3 


82/64 






f 




0-09±g;gi 


f 




S2+°' 01 
U -J z _0.01 


f 




15+ ' 01 






f 




3 4+ 01 

°-^-0.2 




















Q 




13+.° 


Q 




i-5±8:? 




















P20181.07 


30 


A 




14+ 000 
u - i ^-o.oi 


R 




05+°'°° 


R 




o.08i°;°; 


R -- 


- u.io_ 00 


A 




6 3+°' 2 
D.O_o 2 


X 


10- 


3 


97/61 






f 




10+° 01 
u - ±u -0.00 


f 




22+°' 01 
u - zz -0.00 


f 




n 00+0.01 

U.38_ 00 


f = 


- 47+ 01 
-0.01 


f 




7 +0.1 
°-'-0.2 




















Q 




8.6+ 3 2 j 


Q 




13+2 


Q-- 


- 1 R+ ' 1 
















P20181.08 


20 


A 




17+ - 01 


R 




07+ ' 01 
u - u -0.01 


R 




15+ 001 

u - io -o.oi 






A 




5 9+ ' 1 


X 


10- 


3 


64/59 






f 




n+ 001 


f 




3fi+ 002 


f 




u - 4± -o.oi 






f 




0+0.2 

d -°-0.3 




















Q 




4 9+1.6 


Q 




1 0+0.4 

L - 6 -0A 
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Table A5. Fits to the 10-20 keV PSD of GX 339-4 (1999 hard state). 



Obs. 


/max (Hz) 


zfc-Lori 


QPOi 


QP0 2 


QPO3 


zfc-Lor2 


X 2 /DoF 


P40108.03 


8 


* = 0.021°;°° 

- or + 0.05 

/ = 0.35+^ ^ 


i? = 0.l6l°;°i 

„ +0.30 

/ = 0.961^ 

Q = o.7l°; 4 2 








53/38 


P40108.04 


15 


A -ft ii+O.OO 

— u.ii_ 01 

p r. i ^ + 01 


R — n 14+001 

p r. rr + D 02 

/ = 0.55!£^ 
Q = 1.8l°i 






A = 3.9l°'| x 10" 3 

/ = 4.0!^ 


110/95 


P40108.05 


25 


A = 0.20+°°; 

v u — 0.01 


R = 0.07+°°° 

c u ' —0.01 


R = 0.09+?:™ 

1 u ^ — 0.01 


R = 0.05+°°° 


A = 9.6+°^ x 10" 3 


56/54 






/ = 0.12±g;gl 


/ = 0.28l°;°? 
= 111* 


/ = 0.40l°;° 
Q = 6.ll^ 


/ = 0.64l°;° 2 2 
Q = lll 9 3 


/ = 2.7l°;i 




P40108.06 


8 


A = 0.42±°;°* 

/ = 0.051°;°; 


fl = 0.181°;°^ 

/ = o.i4l°;°? 
Q = i.ol°? 






A = 10.7+°'° x 10- 3 

/ = 2.5l°; 2 


53/36 


P40108.07 


5 


A = 0.431°;°;? 

/ = 0.051°;°? 








A = 8.4+^ x 10- 3 

/ = i.5±S;S 


20/19 
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